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CHAPTER  1 


INTRODUCTION 


1.1.  Summary  of  Previous  Results 


i 


In  two  previous  papers  ([5]  and  [6]),  stochastic  capacity  ex- 
pansion Models  I and  II  were  introduced.  Both  models  permit  two  types 
of  facilities:  permanent  and  temporary.  Permanent  facilities  represent 
the  means  by  which  demand  is  normally  accommodated,  while  temporary 
facilities  represent  the  extraordinary  measures  taken  in  order  to 
accommodate  excess  demand  prior  to  an  expansion  of  permanent  facilities. 
Examples  of  temporary  facility  usage  include  cases  of  backlogging,  over- 
loading and  jobletting. 

Models  I and  II  were  constructed  under  the  following  assumptions: 

a.  The  demand  process  and  costs  are  time-stationary. 

b.  The  demand  process  can  be  characterized  as  the  difference  of  two 

k 

simple  Poisson  processes  (without  truncation) . 

c.  All  costs  are  independent  of  the  amount  of  permanent  capacity 
available. 

d.  Permanent  facility  operating  costs  are  proportional  to  the  amount  of 

kk 

permanent  facility  capacity  actually  utilized. 


The  untruncated  demand  assumption  continues  to  permit  departure  when  the 
demand  level  is  zero.  The  viability  of  this  assumption  is  discussed  at 
length  in  Section  1.7  of  [5]. 

It  is  demonstrated  in  Sections  1.5  of  both  [5]  and  [6]  that  permanent 
facility  operating  cost  terms  which  are  proportional  to  the  permanent 
capacity  available  are  also  permissible. 
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Model  I,  introduced  in  [5],  addresses  the  case  where  temporary  facilities 
are  non-mo dular  with  usage  costs  depending  solely  on  the  level  of  excess 
demand  being  served.  Model  II,  introduced  in  [6],  treats  the  case  of 
modular  temporary  facilities  which  are  typically  available  in  a fixed 
increment  size  (a  module)  and  incur  instantaneous  installation  and  removal 
charges  in  addition  to  normal  usage  costs.  The  objective  in  each  model 
is  the  minimization  of  expected  total  discounted  costs.  The  expansion 
decision  strategies  considered  are  of  the  recurrent  form  (X,K) , where 
K denotes  a limit  on  temporary  facility  usage  and  X+l  denotes  the 
permanent  facility  expansion  size.  That  is,  whenever  excess  demand  equals 
K,  a permanent  expansion  of  size  X+l  is  implemented  upon  the  arrival  of 
an  additional  demand  unit. 

In  each  model,  the  assumptions  (a)  - (d)  allowed  the  derivation 
of  a single  finite  linear  system  of  expected  cost  equations  for  fixed 
values  of  X and  K.  In  the  case  of  Model  II,  the  derivation  of  the 
linear  equation  system  was  preceeded  by  a demonstration  that  an  optimal 
module-removal  policy  is  of  the  form:  "remove  a module  only  if  at  least 

s*  units  of  unused  capacity  will  remain  available".  An  algorithm  was 
given  to  determine  the  optimal  parameter  s*. 

By  parameterizing  the  expected  cost  equation  systems  in  the  ex- 
pansion variable  X,  it  was  shown  that  linear  equation  solution  techniques 
can  be  used  to  derive  expected  cost  functionals  Cq(*,K)  for  each  fixed 
value  K.  It  was  also  shown  that  an  optimal  expansion  size  X*(k)  + 1 
satisfies 


Cq(X*(K),  K)  = min{CQ(X,K)  : integer  X > K>  . 


2 


Algorithms  were  given  in  [5]  and  [6]  for  recursively  computing  these 
functionals  over  all  feasible  values  for  K.  The  form  of  the  func- 
tionals was  identical  for  both  Models  I and  II. 

This  paper  consists  of  two  self-contained  parts.  In  Chapter  2, 
the  determination  of  optimal  expansion  sizes  for  recurrent  Models  I 
and  II  is  considered.  In  Chapter  3,  a revised  non-recurrent  model, 
permitting  the  relaxation  of  assumptions  (b)  - (d)  above,  is  intro- 
duced . 


1.2.  Notation 

-S  ( t ) will  denote  the  demand  at  time  t >_  0.  Given  Poisson 
arrival  rate  A^  > 0 and  Poisson  departure  rate  A^  > 0,  A = A^  + A^ 
will  denote  the  combined  Poisson  event  rate , with  arrival  probability 
p = A^/A  > 0 and  departure  probability  q = A2/A  > 0(p+q  = 1)  [3]. 
is,  whenever  an  event  occurs,  the  event  is  a unit  demand  increase  with 
probability  p and  the  event  is  a unit  demand  decrease  with  probability 

q* 

As  introduced  earlier  in  this  chapter,  K will  benote  the 
temporary  facilities  usage  limit  and  X+l  will  denote  the  permanent 
facilities  expansion  size,  g will  denote  the  expansion  cost  function, 
where  g(X)  represents  the  cost  of  an  expansion  of  size  X+l.  r will 
denote  the  continuous  discount  rate;  0 < r < 1. 

k will  denote  the  spares  level:  the  permanent  capacity  available 
less  demand.  Thus,  whenever  k _>  0,  the  permanent  facilities  available 
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CHAPTER  2 


COMPUTING  OPTIMAL  EXPANSION  SIZES  FOP  THE  RECURRENT  MODELS 
2.1.  Introduction 

[5]  and  [6]  demonstrate  that  for  a fixed  temporary  facility 
usage  limit  K 6 K (where  k denotes  a nonnegative  integer  set  of 
feasible  limits  with  0 6 k) , an  optimal  permanent  expansion  size 
X*(K)  + 1 satisfies 


where 


C0(X*(K),  K)  = min{C0(X,K)  : integer  X > K}  (2.0) 


C0(X,K)  = ^7l~+  ^ ZK  , X > K ; (2.1) 


a 

) is  nonnegative  and  nondecreasing  over  «;  (2.2) 


0 < Z 


— ~ < 1,  given  a = Aq(A+r)_1 


- 1 2 

and  (3  = Ap(A+r)  (Z  satisfies  OZ  - Z + (3  = 0) ; (2.3) 


Ck(X,K)  = l C0(X,K)  , where  C^X^K)  denotes  the 
expected  total  discounted  costs  starting  from  initial  spares 
level  k > 0 (with  no  temporary  facilities  initially  in  use) 


Algorithms  for  computing  { 0 (K ) : K £ «}  are  given  in  [5]  and  [6]. 


5 





r 


. 

. 


] 

i 


j 


and  using  permanent  expansion  size  X+l  (X  > K)  when  the 
temporary  facilities  usage  limit  has  value  K € k;  (2.4) 


Z 


k 


where  denotes  the  first  time  that  the 

net  demand  increases  by  an  amount  k > 1.  (2 


This  chapter  treats  the  above  problem  under  two  basic  assumptions: 
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Assumption  2.1.  g is  nonnegative  and  nondecreasing  over  [o.co). 

Assumption  2.2.  There  exists  a known,  finite,  nondecreasing  function 
U( • ) over  k such  that  X*(K)  € [K,  U(K)],  K £ k. 

Assumption  2 . 1 is  self-explanatory  and  anything  to  the  contrary 
would  be  unrealistic.  Assumption  2.2  is  used  to  insure  that  the  search 
for  X*(K)  can  be  restricted  to  the  finite  interval  [K,  U(K)].  In 
practical  situations,  unlimited  expansion  is  obviously  infeasible.  The 
following  lemma  provides  mathematical  justification  for  this  point  of 
view. 


I 

i 


i 

j 


! 

' . 


Lemma  2. 1.  If  there  exists  A > 0 such  that  g(X+A)  - g(X)  > ? > 0 
for  all  X > 0,  then  Assumption  2.2  is  satisfied  with 


U(K)  = max{  K , U’(K),  [|  (Z-"..^(K) 


,A 


K ^ ^ - 0(K))  + 1]}  , K € < 
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where  U'(0)  = 0,  and 

U'(K)  = max { U( K 1 ) : K ' € x ~ [K,w)}  , K £ * ~ {0} 


( the  brackets  denote  integer  truncation) . 


Proof.  The  inclusion  of  U'(»)  insures  that  U(*)  is  nondecreasing 
over  k.  Suppose  X>U(K).  Then 


+ g(K) ) 

1 r,K 


- 0(K))  + 1 


=»g(X)  > g(X)  - g(0)  > |[|]  > 2 + ffW)  . 0(K) 

z - z 

=M(K)  + g(X)  (P(K)  * S(K)) 

z'1  - zK 

^(-K-l  + gffl  -^(K)  g(K) 

z-1  - zx  z'1  - zK 

^C0(X,K)  >C0(K,K)  . 


□ 


The  above  leirana  essentially  states  that  unlimited  expansion  is 
never  optimal  if  the  marginal  costs  of  expansion  are  bounded  away  from 
zero. 

Denote 


C(X,K) 


&K)  + 8(X) 
z-1  - zx  ’ 


x > 0,  K 6 k . (2.6) 
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Also,  by  (2.0)  and  Assumption  2.2,  X*(K)  is  equivalently  determined  by 


C(X*(K),  K)  = min{C(X,K)  : integer  X £ [K,  U( K ) ] ) . 


(2.8) 


Lemma  2.2.  For  X > 0f  C(XJ*)  is  nondecreasing  over  x. 


Proof.  Let  K,K'  £ < with  K >K'.  By  (2.2),  0(K)  >P(K').  Hence^ 


C(X,K)  = 


z-1  - zx  - z-1  - zx 


= C(X,K') 


Definition  2 . 1.  The  function  g is  convex  over  the  interval  [a^b]  if 
g(nX  + (l-n)X')  < (_ig( X)  + (1-p)  g(X')  for  all  real  X,X’  £ [a,b], 

0 < p.  < 1.  The  function  g is  concave  over  the  interval  [a^b]  if 
-g  is  convex  over  the  interval  [a^b].  The  function  g is  discrete 
convex  over  the  interval  [a^b]  if  g(uX  + (l-u)X')  < ug(X)  + ( 1-p)  g(X') 
for  all  integers  y.X1  £ [a^bj^  0 < p < 1.  The  function  g is  discrete 
concave  over  the  interval  [a?b]  if  -g  is  discrete  convex  over  [a^b]. 
Recall  that  g(0)  represents  the  cost  of  a permanent  expansion 
of  size  1.  Hence?  a restriction  of  g to  convexity  over  [0,oo)  does 
not  preclude  the  existence  of  a fixed  (or  "setup")  charge  in  the 
expansion  cost  function.  Notice  that  if  g is  convex  (concave)^  then 
g is  discrete  convex  (discrete  concave). 
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In  this  section,  the  determination  of  an  optimal  permanent  expan- 


sion size  X*(K)  + 1,  for  a fixed  value  K £ K}  is  considered  for  arbitrary 
nonnegative,  nondecreasing  functions  g.  The  approach  used  is  based 
on  Howard's  well-known  Policy  Improvement  Technique  of  dynamic  programming 
[ l].  In  the  context  of  the  problem  at  hand,  the  Policy  Improvement  result 
can  be  paraphrased  as  follows  for  fixed  K £ k: 

"Using  expansion  size  X+l  is  more  costly  than  using 
expansion  size  X'+l  if  and  only  if,  using  X+l  is 
more  costly  than  using  X'+l  for  the  first  expansion  and 
X+l  thereafter." 

Theorem  2.5  is  the  precise  statement  of  this  result. 

Definition  2 .2  . Let  ^k(K),  k > 0,  denote  the  expected  discounted 
operating  costs  until  the  first  permanent  expansion,  when  the  initial 
spares  level  is  k > 0 (and  no  temporary  facilities  are  initially  in 
use)  and  the  temporary  facilities  usage  limit  is  K £ K. 


Lemma  2.3.  For  K £ k. 


ck(X,K)  = Gk(K)  + (g(X)  + Cx_k(X,K))  ZK+k+1  , k > 0,  X > K 


Proof.  Starting  from  spares  level  k,  the  first  expansion  occurs  whenever 
the  net  demand  first  increases  by  an  amount  K+k+1;  this  time  is  given 
by  i anc*  zK+*c+^  _ E[e  K+k+1  j . xhe  expecte(j  discounted  costs 


over  the  time  interval  [0,  TR+k+1)  are  Gk(K).  At  time  0^  an 

expansion  is  undertaken  at  cost  g(X)  and  the  expected  discounted  costs 
thereafter  are  GXK(X,K),  relative  to  . Hence, 

ck(x,K)  = GJK)  + E[(g(X)  + Cx_k(X,K))  e'rTK+k+1] 

= Gk(K)  + (g(X)  + Cx_k(X,K))  ZK+k+1  . □ 

Lemma  2.4.  For  K ■ E K, 

„K+k+ 1 

Ck(X,K)  = Gk(K)  + — ^ (g(X)  + Gx_k(K))  , k > 0,  X > K . 

I"  o 

Proof.  From  the  previous  lemma, 

<WX>K>  = ^X-k(X>K)  + (&(X)  + CX_X(X>K))  Z 

Iterating  the  above  equation  yields 

Cx_k(X,K)  = Gx_k(X,K)  + (g(X)  + gx  k(x,k))  z^x+1) 

i=  1 

+ (g(X)  + CXk(X,K))  Zfi(X+1),  l = 1,  2,  ... 
Letting  },  ->  oo  in  the  above  equation  gives 

■ CX-K(X'K>  = + (8(X)  + Gx.k(x,K))  z Z^X+1)  . 

I 

* 

I 
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Substituting  this  expression  into  that  of  the  previous  lemma  yields 


ck(X,K)  = Gk(K)  + (g(X)  + Gx_k(K))  (1  + £ z1(X+l))  ZK+k+1 


= Gk(K)  + ZK+k+1(g(X)  + (K))  Z (Z^+V 

i=0 

zK+k+l 

= VK)  (g(x)  + Gx_k(K))  . 


Theorem  2.  5.  (Policy  Improvement).  For  X,X'  >K  £ k 

C(X',K)  |.  |c(X,K)  g(X')  ♦ Cx,_k(X,K)  | =|  g(x)  + CX  K(X,X)  . 

Proof.  Wo  shall  prove  the  the Or™  using  the  < Inequality;  the  other 
two  cases  can  be  proved  in  an  analogous  manner. 

We  shall  first  prove  sufficiency  (<=) . Let  h = ZX'+1.  By 
Lemma  2.3  and  the  hypothesis 


cx’-k^X'K)  = WK)  + (s(x)  + cx.K(xpK))h 

> °X --K(K)  + (s(X')  + C (X,K))h 


(2.9) 


Iterating  (2.9)  £+1  times  yields 


CX'-K(X>K)  >GX'_k(K)  + (g(x>)  + G , (X/))  Z h1 


+ (8^X')  + cx'-k^k))  hU\  * = 0,  1,  ... 


(2 . 10) 


Since  (2.g)  is  applied  at  each  iteration,  the  righthand  side  of  (2.1o) 
decreases  with  £.  Hence,  strict  inequality  is  maintained  in  the  limit, 
yielding  as  i -» » that 


i 


V-K<M>  > GX'-K<K>  +1Th  <6<X')  * °X'-K<K»  ' CX'-K<X’>K>  - 

where  the  equality  follows  from  Lemma  2.4.  By  (2.7),  the  above  inequality 
implies  that  C(X,K)  >C(X',K). 

The  proof  of  necessity  (=>)  is  by  contradiction.  Suppose  that 
C(X',K)  < C(X,K),  but  contrary  to  the  theorem, 

g(x')  + Cx,_k(X,K)  > g(X)  + Cx_k(X,K)  . (2.11) 

X '+1 

Let  h = Z . By  Lemma  2.3  and  the  above  inequality^ 

WM)  = Gx'.rCK)  + (s(X)  + Cx_K(X'K))h 

< Gx,_k(K)  + (g(X')  + Cx,_K(X,K))h  . 

Iterating  this  expression  to  the  limit  gives 


cx._k(x,k)  <gx,_k(k)  + (s(x')  + <^,.K(K))V(i-h)  - cx,_k(x*,k)  , 

\ 

1 


which  by  ( 2.7),  implies  that  C(X,K)  < C(X',K).  Since  this  contradicts 
the  hypothesis,  we  must  conclude  that  the  supposition  (2.1$  is  false, 
which  completes  the  proof.  □ 
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Using  (2.7) ^ the  Policy  Improvement  Theorem  can  be  restated  as 


follows. 


Corollary  2.6.  For  X^X'  > K £ k, 


C(X',K)  = C(X,K)  «=»g(X')  - g(X)|=  C(X,K)  (Z^Z*') 


Proof. 


(<)  (<) 

C(X',K)  j=  jC(X,K)  «=*g(X')  - g(X)  j = jCX  K(X,K)  - CX,_K(X,K) 


= C(X,K)  (Z***  - ZK^'-K) 


= C(X,K)  (ZX  - ZX') 


Corollary  2.7.  Let  X > K £ k.  For  0 < x < U(K)-X, 


C(X+x,K)  j=jc(X,K)*=>  g( X+x)  - g(X)  j = ) C(X,K)  ZX(  1-ZX)  (2.12) 


For  0 < x < X-K, 


C(X-x,K)  J =jc(X,K)  <=*  g(X)  - g(X-x)  |=  |c(X,K)  ZX(Z_X-1)  . (2.13) 
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Proof.  (2.12)  follows  from  Corollary  2.6  upon  setting  X*  = X+x  and 
(2.13)  follows  likewise  upon  setting  X'  = X-x.  □ 

Using  Corollary  2.7,  the  Policy  Improvement  Algorithm  for  finding 
X*(K)  can  be  stated  in  two  different  ways.  Initializing  the  algorithm 
at  the  beginning  of  the  interval  [K,  U(K)]  and  using  (2.12)  gives: 

Algorithm  C (Policy  Improvement). 

(1)  XQ  = K;  i = 0. 

(2)  Compute  C(X^,K). 

(3)  Find  x*  = min{x  : g(X.+x)  - g(X.)  < C(X£,K)  Z L(l-Zx),  0 < x < U(K)-X  ). 
If  x'  does  not  exist,  then  X*(K)  = X^  (stop). 

(4)  X_.^  = X^+xj  i «- i+l;  go  to  step  (2). 

The  sequence  {X^j  produced  by  Algorithm  C^  is  increasing  and  (C(Xi,K) ) 
is  decreasing.  Since  each  (integer)  X.  £ (K,  U(K) ],  the  algorithm  is 
finite. 

Initializing  the  algorithm  at  the  end  of  the  interval  [K,  U(K) j 
and  using  (2.13)  gives: 

Algorithm  C^  (Policy  Improvement). 

(1)  XQ  = U(K) } i = 0. 

(2)  Compute  C(X^,K). 

(3)  Find  x'  = min{x  : g(Xt)  - g(X.-x)  > C(X1,K)  Z 1(Z'X-1),  (Xx^Xj-K). 

If  x'  does  not  exist,  then  X*(K)  = X^  (stop). 
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(4)  X = X^-x i *- i+1;  go  to  step  (2). 

The  sequence  {X^}  produced  by  Algorithm  C2  is  decreasing,  as  is 
{C(Xi,K)}.  Since  each  (integer)  X^  £ [K,  U(K) ],  the  algorithm  is  finite. 

It  should  be  noted  that  Corollary  2.6  can  also  be  used  to  state 
an  algorithm  that  permits  an  arbitrary  initial  choice  XQ  £ [K,  U(K)j. 

In  this  case,  the  sequence  (C(X^,K) } will  again  be  decreasing,  and 
hence  the  algorithm  is  finite  since  (X^)  will  be  a nonrepeating  sequence 
on  a finite  set.  However,  the  sequence  (X^}  will,  in  general,  not  be 
monotone;  this  is  often  impractical.  For  this  reason,  Algorithms  and 
would  be  preferred  in  practice. 

While  Policy  Improvement  provides  finite  algorithms,  it  is,  in 
the  worst  case,  no  better  than  an  exhaustive  enumeration  of  the  interval 
[K,  U(K)].  However,  the  Policy  Improvement  Theorem  can  be  used  to  obtain 
conditions  on  g which  allow  more  efficient  search  techniques.  The 
next  section  addresses  this  subject. 

2.3.  Computing  X*(K) : Unimodal  Properties.  Interval  Bisection  and 

Monotonicity  for  Discrete  g. 

In  this  section,  the  Policy  Improvement  Theorem  is  used  to  derive 
conditions  on  g which  abet  more  efficient  search  techniques.  General 
conditions  are  derived  for  g which  will  guarantee  that  C(.,K)  is 
unimodal  for  fixed  K £ k. 
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Definition  2.3.  For  fixed  K £ Ky  is  unimoda 1 over  the  interval 

U(K)]  if  there  exists  X'  £ [K,  U(K)]  such  that: 

(i)  C(X,K)  > C(X+1,  K),  K < X < X';  and 
(ii)  C(X,K)  < C(X+1,  K),  X*  < X < U(K). 

Definition  2.4.  For  K £ Kf  let  X^(K)  and  x,2(K)  denote  the  sets 
given  by: 

Xl(K)  = {X  : g(X+l)  - g(X)  > C^K)  ZX(1-Z),  K < X < U(K)  } , 


X2(K)  = {X  : g(X)  - g(X-l)  < C(X,K)  ZX(Z-1-l)j  K < X < U(K)J 


Let  X^(K)  and  Xg(K)  be  given  by 


XX(K)  = 


x2(x)  = 


m. 


if  Xl(K)  = 0 


min(X  : X £ Xl(K)),  if  X]_(K)  ^ 0 

K,  if  Xg(K)  = 0 

tnax(X  : X £ xg(K)),  if  ^(K)  ± 0 


The  following  lemma  provides  alternative  definitions  for  the  above  two 


sets. 
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Lemma  2.8.  For  K € K and  X = [K,  U(K)], 


X 6 X (K)  C(X,K)  < C(X+1,  K) 


X € X2(X)  «=*C(X-1,  K)  > C(X,K) 


Proof.  Follows  from  Definition  2.4  and  Corollary  2.7. 


Lemma  2.9.  For  K € K. 


K < X < XX(K)  -=>  C(X,K)  > C(X+1,  K)  =>  C(X,K)  > C(Xl(K)  , K)  , 


Xg(K)  <X  < U(K)  -*C(X-1,K)  < C(X,K)  =»C(X2(K),  K)  < C(X,K) 


Proof.  From  the  previous  lemma  and  Definition  2.4, 


K < X < Xl(K)  =»  X i x^K)  <=?>  C(X,K)  > C(X+1,K) 

X2(K)  < X < U(K)  =>X  4 ^(K)  «=»  C(X-1,  K)  < C(X,K) 


Thus,  if  X^(K)  / j6,  then  X^(K)  is  the  first  point  of  increase 
for  C(  • ,K)  over  [K,  0(K)  ] (i.e.,  C(X^(K)  + 1,  K)  >C(X1(K),  K)). 

Similarly,  if  X^K)  / j6,  then  Xg(K)  is  the  last  point  of  decrease 
for  C( • ,K)  over  [K,  U(K) ] (i.e.,  C(X2(K)-1,  K)  > C(X2(K),  K)). 
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Theorem  2.10  (Monotoniei ty) 


are  each  nondecreasing 


over  k 


Proof.  Let  K,K ' € k with  K < K 


then  tne  result 


follows  trivially.  Therefore,  assume  ttiat  X^K ')  < U(K)-1 
U(K')-1  =3>  X (K1)  ^ 0.  Hence  using  Definition  2.4  and  Lemma  2.2 


which  implies  that  X^K')  € X^K),  since  X^K')  € [K ',  U(K)-l]  [K,U(K)-1] 

Therefore,  X^K)  ^ 0 and  X^K')  > min{X  : X eX^K)  ) -=  X (K). 

If  X2(K)  < K',  then  the  result  for  X2(  • ) follows  trivially. 
Therefore,  assume  that  Xg(K)  > K'  + l =»  X2(K)  > K+l  =*  X_(K)  f 0.  Hence, 
using  Definition  2.4  and  Lemma  2.2 


8(X2(K))  - g(X2(K)  - 1)  <C(X2(K),  K)  Z 


< C(Xg(K),  K ')  Z (Z  -1) 


which  implies  that  X (K)  € X^K'),  since 
= [K'+l,  U(K')  ].  Therefore,  )C(K')  ^ 0 


X2(K)  6 [K'+l,  U(K)  1 
and  X2(K)  < max{X:X  £ X^K')} 


The  above  theorem  is  independent  of  the  functional  form  for  g 
andj  in  that  context^  it  is  worthy  of  restatement.  The  Monotonicity 
Theorem  states  that  within  k:  (a)  points  of  first  increase  for  C(*^K) 

over  [Kj  U(K)]  are  nondecreasing  in  K;  and  (b)  points  of  last  decrease 
for  C(«jK)  over  [K,  U(K)]  are  also  nondecreasing  in  K. 

The  next  series  of  propositions  derive  conditions  on  g that 
lead  to  general  unimoclal  results. 

Definition  2.5.  Let  K £ k.  g will  be  said  to  satisfy  property  if 

g(X+l)  - g(X)  > Z(g(X)  - g(X-l))  , K+l  < X < U(K) - 1 

g will  be  said  to  satisfy  property  A if  g satisfies  property 
for  all  K G k. 

The  above  definition  is  dependent  on  the  value  Zf  as  determined 
by  the  probability  parameters  (A, p,q)  in  (2.3).  However^  the  following 
result  is  independent  of  these  parameters. 

Lemma  2. 11.  Let  K £ k.  If  g is  discrete  convex  over  [K^  U(K)], 
then  g satisfies  property  A^.  If  g is  discrete  convex  over  tO^oo)^ 
then  g satisfies  property  A. 
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Proof,  g discrete  convex  over  [K^  U(K) J implies  that 


g(X+l)  - g(x)  > g( X)  - g(X.l)  > Z(g(X)  - g(X-l)),  K+l  < X < U(K)-1  , 


since  0 < Z < 1.  If  g is  discrete  convex  over  then  g is 

discrete  convex  over  [K,  U(K)]^  for  all  K £ k.  □ 


The  next  two  lemmas  list  equivalent  definitions  for  the  property 


V 


Lemma  2,12:  Let  K £ k.  g satisfies  property  if  and  only  if 


g(X+x)  - g(X)  >jf| 


(g(X+l)  - g(X) ) 


y 


1 < x < U(K)-X,  K < X < U( K) - 1 


(2.14) 


Proof.  To  prove  necessity  (=>),  suppose  g satisfies  property  A^.  Then^ 


g(X+i)  - g( X+i- 1)  > Z( g(X+i- 1)  - g(X i-i-2) ) 


> . . . 

> z1'  L(  g(  X+l)  - g(  x) ) . 


Hence^ 


g(X+x)  - g(X)  = Z (g(X+i)  - g(X+i-l)) 
i=l 


> (g(x+i)  - g(x))  z z1 

i=l 


= (g(x+i)  - g(x))  i ~ . 


To  prove  sufficiency  («=) , let  x = 2 in  (2.14).  Then, 


g(X+2)  - g(X)  >ff-  (g(X+l)  - g(X) ) 


= (1+Z)  (g(X+l)  - g(X)) 


g(X+2)  - g(X+l)  > Z(g(X+l)  - g(x)) 


Lemma  2.13.  Let  K £ k.  g satisfies  property  if  and  only  if 


Z"x-1 


g(X)  - g(X-x)  <~~  (g(X)  - g(X-l))  , 
Z -1 


i < x < X-K,  K+l  < X < U(K)  . (2 .15) 


Proof.  To  prove  necessity  (=>),  suppose  g satisfies  property  y^. 
g(X-i)  - g(X-i-l)  < Z“1(g(X-i+l)  . g(X-i) ) 


C(X,K)  > C(X+1,  K)  , 

C(X,K)  = C(X+1,  K)  , 
C(X,K)  < C(X+1,  K)  , 


K < X < X2(K) 

Xg(K)  < X < Xt(K) 
XX(K)  < X < U(K) 


□ 
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Theorem  2.15 

(Unimodality).  Let 

K £ K,  If  g satisfies  property  A, 

l\.  * 

(i) 

C(-,K) 

is  unimodal  over 

[K,  U(K)1, 

(ii) 

C(X,K) 

= C(X*(K) , K),  X £ 

[X2(K),  Xx(K)],  and 

(iii) 

C(X,K) 

> C(X*(K),  K),  X € 

[K,  U(K)  ] ~ (X2(K),  X^K)]. 

Proof. 

By  Definition  2.4 

g(Xx(K)  + 1)  - g(X1(K))  >C(X1(K),  K) 


zXl(K) 


(1-Z) 


y 


\(K)  / 0 . 


Combining  this  with  inequality  (2.14)  gives 


g(X;(K)  + x)  - g(XL(K))  >C(X1(K),  K)  ZXl;K\l_ZX)  , 


1 < x < u(k)-xl(k),  Xi(k)  / p 


Hence,  by  Corollary  2.7? 


C(X,K)  >C(X1(K),K),  X^K)  < X < U(K),  X^K)  / p 


( 2.16) 


Similar ly,  by  Definition  2.4, 


gCX^K))  - g(X2(K)-l)  <C(X2(K)>K)  ZX2lK)(Z_1-l)  , X2(K)  / 0 . 

Combining  this  with  inequality  (2.15)  gives 


then 


23 


xo(K) 

g(X2(K))  - g(Xg(K)  - x)  < C(X2(K),  K)  Z ^ (Z'X-1)  , 

1 < x <X2(K)-K,  Xg(K)  0 0 . 

Hence,  by  Corollary  2.7, 

C(X2(K),K)  < C(X,K),  K < X < Xg(K),  Xg(K)  0 0 . (2.17) 

We  now  consider  four  possible  cases. 

Case  1:  X^K)  = 0 and  x2(x)  = 0*  Then  X2(K)  = K < U(K)  = XL(X)  and 

(i)-(ii)  follow  from  the  previous  lemma.  ( iii)  is  null  for  this  case. 

Case  2:  Xl(K)  0 0 and  ^(K)  = 0.  Then  Xg(K)  = K < X^K)  and  (i)-(ii) 

follow  from  the  previous  lemma,  (iii)  follows  from  (2.16). 

Case  3-'  X^K)  = 0 and  X^K)  0 0.  Then  Xp(K)  < U(K)  = X^K)  and 

(i)-(ii)  follow  from  the  previous  lemma,  (iii)  follows  from  (2.17). 

Case  4:  Xj_(K)  0 0 and  X2(K)  0 0.  Suppose  X^K)  <Xg(K).  Then  (2.16) 

and  (2.17)  together  imply  that  C(Xg(K),  K)  >C(X1(K),  K)  >C(Xg(K),  K) , 
a contradiction.  Hence,  X^(K)  > Xg(K)  and  (i)-(ii)  follow  from  the 
previous  lemma.  (iii)  follows  from  (2.16)  and  (2.17).  □ 
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By  (ii)  and  ( iii)  of  the  above  theorem,  X^K)  and  Xg(K)  are 
respectively  the  greatest  and  least  elements  in  the  minimizing  set  for 
C(*,K)  over  [K,  U(K)].  Hence,  whenever  g satisfies  property 
X2(K)  + 1 is  the  least  optimal  expansion  size  and  X^Kj+l  is  the  greatest 
optimal  expansion  size.  Combining  this  observation  with  the  Monotonicity 
Theorem  gives 

Corollary  2,16  (Monotonicity).  If  g satisfies  property  A,  then 
(i)-(iii)  of  Theorem  2.15  hold  for  all  K 6 k,  and  the  least  and  greatest 
optimal  expansion  sizes  (1+X^(‘)  and  1+X1(.),  respectively)  are  each 
nondecreasing  over  k. 

Proof.  The  result  follows  directly  from  Theorems  2.10  and  2.15.  □ 

Corollary  2,17.  Let  K £ (cj  if  g is  discrete  convex  over  [K,  U(K)  ], 
then  (i)-(iii)  of  Theorem  2.15  follow.  If  g is  discrete  convex  over 
[0,«),  then  (i)-(iii)  of  Theorem  2.15  hold  for  all  K £ k and  the 
least  and  greatest  optimal  expansion  sizes  (1+Xg(*)  and  1+X^(.), 
respectively)  are  each  nondecreasing  over  k. 

Proof.  The  result  follows  directly  from  Theorem  2.15,  Corollary  2.16, 
and  Lemma  2.11.  □ 
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The  above  propositions  detail  conditions  on  g which  guarantee 
that  C(*,K)  is  unimodal  over  [K,  U(K)].  In  these  cases,  X*(K)  can 
be  found  by  using  Interval  Bisection  ("Binary  Search"): 


Algorithm  C,  (Interval  Bisection,  given  , 


(1) 

x3 

II 

ii 

c3 

s 

• 

(2) 

if 

Xi  " Xu>  then  X*(R)  = 

*£  = Xu  (8toP)- 

(3) 

x = 

: [(X  +X  )/2], 

L U 

(M 

if 

g(x+l) 

- g(X)  > C(X,K) 

x 

Z.  ( 1-Z)  , then  reset 

reset  X = 
u 

1 x+1. 

(5) 

Go 

to  step 

(2). 

= X;  otherwise, 


Using  Algorithm  C^,  the  search  interval  is  essentially  halved  at  each 
iteration  and  no  more  than  1 + log^UfK)  - K)  iterations  are  required. 
Note  that  by  Corollary  2.7,  the  "if"  condition  of  step  (4)  is  equivalent 
to:  "If  C(X,K)  < C(X+1,  K) ",  Thus,  the  iterate  X is  only  reset  to 

Ju 

X if  C(X',K)  >C(X'+1,K)  for  X'  < X-l.  Hence,  the  algorithm  computes 
X2(K)>  the  least  optimal  expansion  size  (less  unity).  By  removing  the 
equality  portion  of  the  "if"  condition  in  step  (4),  the  algorithm  will 
alternatively  compute  X^(K),  the  greatest  optimal  expansion  size  (less 
unity).  Obviously,  the  algorithm  can  be  easily  modified  to  compute  both 
X^(K)  and  Xg(K),  thereby  giving  the  entire  set  of  optimal  expansion 
sizes  for  K : [Xg(K)+l,  X^Kj+l]. 

In  view  of  Corollary  2.16,  some  additional  efficiency  can  be 
realized  in  Algorithm  C^  whenever  it  is  used  to  recursively  compute  X*(») 
over  k.  To  illustrate,  suppose  that  X'  and  K are  consecutive  elements 
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of  < with  k'  < K.  Then  X*(K ')  is  a lower  bound  on  X*(K)  and 
additional  efficiency  inay  be  realized  by  alternatively  initializing  the 
lower  iterate  X,  = max{K,  X*(K')}  in  order  to  compute  X*(K). 

The  final  theorem  of  this  section  provides  a general  result  that 
may  be  used  to  construct  more  efficient  search  algorithms  for  a number 
of  other  cases  where  Property  cannot  be  ascertained.  It  demonstrates 
that  further  restriction  of  the  searcti  interval  may  be  possible,  given 
known  bounds  on  the  marginal  costs  (e.g.,  g a step  function).  This 
result  can  be  used  for  arbitrary  g is  conjunction  with  Algorithm 
of  the  previous  section. 


Theorem  2.18.  Let  K £ K.  Suppose  0 < f;  ^ < g(X+l)  - g(X)  < > f°r 

X £ [K,  U(K)].  Then  X*(K)  £ [X^K),  U^K)],  where  U (K)  = min[U(K),  X (K) 

+ [«2  r^i-z)-1]). 


Proof.  By  Lemma  2.9,  it  sufficies  to  show  that  X > U^(K)  implies 
X /.  X*(K)  whenever  U^'K)  <U(K).  When  U^K)  <U(K), 


X > Ut(K)  =>  X - XX(K)  = x > [t2  |'11(1-Z)'1]  + 1 


g(X)  - g(XL(K))  = g(X1(K)  + x)  - g(X1(K)) 


> 5j_x 


. . 1-Z* 


- f‘2  1-Z 


> -xrf-  (g(X1(K)  + l)  - g(XL(K))  . 
By  Corollary  2.7,  this  implies  that  C(X,K)  >C(X^fK),  K) . 
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2.4.  Computing  X*(K) : Function  Iteration  and  Monotonicity  for 

Concave  g 

In  this  section,  the  determination  of  X*(K)  is  considered  for 
the  important  case  where  g is  concave.  Assumption  2.2,  introduced  in 
Section  2.1,  will  be  retained.  However,  Assumption  2.1  will  be  strengthened: 

Assumption  2.1 g is  positive,  increasing  and  continuously  differ- 
entiable. 

Under  Assumption  2.1',  C(*,K)  is  continuous  for  each  fixed 
value  K > 0.  Unlike  previous  sections  of  this  chapter,  X will  be 
treated  as  a continuous  variable  here.  However,  X*(K)  + 1 will  continue 
to  denote  the  optimal  integer  expansion  size.  Notice  that  the  positivity 
of  g is  not  impractical  since  g(0)  has  been  defined  as  the  cost  for 
a permanent  expansion  of  size  1. 

Definition  2.5. 


f(X,K)  . -i-rj  ln(  1 . In  Z'1)  .1, 


In  Z 


X > 0,  K € K. 


The  following  lemma  demonstrates  that  the  function  f is  well- 
defined  under  Assumption  2.1 


Lemma  2.19.  f(X,K)  > - 1,  X > 0,  K € k. 


28 


29 


Thus,  if  X < +»  is  a global  minima  of  C(*,K)  over  [0,oo), 
then  X must  be  a fixed  point  of  f(*,K);  that  is,  X = f(X,K).  The 
behavior  of  f with  resect  to  K is  given  by  the  following  lemma. 


Lemma  2.21.  f(X,  •)  is  nondecreasing  over  k,  X ■>  0. 


Proof.  Let  K,  K ' £ * with  K'  < K.  Let  X > 0.  By  (2.2)  and  positivity. 


i + ln  z- 1 < i + in  z-1 

+ g'(x)  - * g'(x) 


Taking  logarithms  and  multiplying  both  sides  by  ( ln  Z- ^ > 0 yields 
f (X,K  *)  + 1 < f(X,K)  + 1.  □ 


When  g is  concave,  f is  also  monotone  with  regard  to  X. 


Lemma  2,22.  If  g is  concave,  then  f(*,K)  is  increasing,  K € K. 


Proof.  Let  K € K,  X > 0 and  e > 0.  When  g is  concave,  g'  i« 
nonincreasing.  Hence 
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ggOHl 


s(x)  < 8(x+e)  - gf(x-re)e  and  g'(X)  > g’fXt-e)  >0 


=*> 


< J&ltiL  - s,(x^<.£ijLt±l 

g’(x)  g'(X) 


< s(x+e) 

- g '(X+e) 


^-St-X)  + 01  K) 

g’(x) 


=£  v = 


/ 1 ,.  s(*~)  * m ln  r1' 

g (x+e) 


> 1 


^ ln  v > 0 


=>  f(X+e,  K)  - f(X,K)  = — - — - In  v > 0 . 

ln  Z 


Definition  2.6.  For  K C K,  condition  £^  will  be  said  to  hold  if  g 
is  concave  and  if  there  exists  ek  > 0 such  that 

f(X,K)>X,  X € (0,  tR] 

and 


f(',K)  is  concave  over  [e  a,)  . 


Condition  £ will  be  said  to  hold  if  holds  for  alL  K € «. 

lx 

Later  in  this  section^  it  will  be  demonstrated  that  condition  £ 
holds  for  several  major  classes  of  concave  functions.  The  following 
lemma  lists  several  properties  of  the  above  conditions. 
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Lemma  2.23,  Suppose  g is  concave. 


o 

(i)  If  f(-,K)  is  concave,  then  (l£)  is  satisfied,  K £ K. 

( ii)  If  (Z^,)  is  satisfied,  for  some  K'  £ K,  then  (Z^)  is  satisfied 

for  all  K > K',  K £ K. 

( iii)  If  f( 0,0)  > 0,  then  (Z^)  is  satisfied  for  all  K £ k. 

2 

( iv)  If  (Z^)  is  satisfied  for  all  K £ k and  f(0,0)  > 0,  then  Z 

holds. 

(v)  If  f ( • ,K)  is  concave  for  all  K £ < and  f(0,0)  > 0,  then 

Z holds. 


2 

Proof.  f(-,K)  concave  over  fo,c»)  obviously  satisfies  (Z^),  so  ( i) 
follows  trivially.  By  Lemma  2.21,  f(X,K)  > f(X,K')  for  K > K 'y  K',  K £ «. 
Hence,  f(X,K’)  > X,  X € (0,  €jc , ] =^>  f(X,K)  > X,  X £ (0,  eK , ] , K > K1,  K £ k; 
this  proves  ( ii) . By  continuity,  f(0,0)  > 0 implies  £ > 0 such 

that  f(X,0)  >X  for  X £ [0,  hence  (iii)  follows  from  (ii).  ( iv) 

follows  from  (iii).  (v)  follows  from  ( i)  and  ( iv) . □ 

; 

The  next  lemma  implies  that,  in  typical  applications,  f will  be 
positive  whenever  g is  concave  (and  hence,  (Z^)  will  be  satisfied 
for  all  K £ k --  see  (iii)  above). 

» 

Lemma  2.24.  If  < 1 - ■ | -,  then  f(0,0)  > 0. 

< 

I I 


- 


g(o))/g'(o) 


The  natural  logarithm  expansion 


8(0))/g'(0) 


then  (r  ) is 


Lemma  2.  25.  If 


Proof.  This  result  follows  directly  from  Lemma  2.24  and  Lemma  2.23^(iii) 


In  practicej  when  g is  concave^  it  is  reasonable  to  expect  the 
hypothesis  of  Lemma  2.24  to  hold.  This  will  now  be  demonstrated.  Using 
(2.3)}  Z can  be  parameterized  in  (rf  A,  p) : 


z(r,A,p) 


( A+  r ) - V(A+r)2  - 4P(l-p)A2 

2A(  1-pi 


(2.18) 


Although  the  proof  of  the  following  lemma  is  analytic^  the  result  can 
also  be  obtained  by  probabilistic  arguments  using  (2.5). 


Lemma  2,26.  Z decreases  with  r increases  with  A,  and  increases  with  p. 


Proof.  Denote 


v = (A+r)  - V(  A+r)2  - 4p(  l-p)A2  > 0 


z = %'filp)7  -^vl/2  » (A+r)  - 2Mi-p)z 

=*  v = (A+r)2  - 4A(A+r)(  l-p)Z  + 4A2(l-p)2  Z2  . 


These  relations  will  be  used  extensively. 

Taking  the  derivative  of  Z with  respect  to  r and  multiplying 

. . 1/2  -1/2  . 

by  1 = v v yields 

sgn(^|)  = sgn( v1'  2 - (A+r))  = sgn(-2A( l-p)Z)  , 


so  dZ/dr  < 0. 
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Taking  the  derivative  of  Z with  respect  to  A and  multiplying 

, . 1/2  -1/2  . .. 

by  1 = v v yields 


dZ  2 

sgn(-j^)  = sgn((A+r)"'  - A(  A+r) 

2 

= sgn((A+r)  - A(A+r) 
= sgn(  ( A+r)£  - (A+r)^ 
= sgn(2Ar( 1- p) Z)  f 


1/2. 
rv  ' ) 

r(  (A+r)  - 2A(1-P)Z)) 
2Ar( l-p)Z) 


so  dZ/dA  > 0. 

2 

Recall  (see  (2.5))  that  Z satisfies  G Z -Z+J3  = 0.  Taking  the 
derivative  of  Z with  respect  to  p and  multiplying  by  1 = v1(/ 2 v"  ^ 2 
gives 


sgn(^)  = sgn(  -v  + (A+r)  v1//<:’  + 2A2(  1- p)  ( l-2p) ) 

= sgn(-2A(l-p)Z2  + (A+r)Z  + A(l-2p)) 

= z" + z * ^T21) 

- ^ z2  + 2 - + g-z2)) 

= sgn((-OZ2  + Z - p)  + P(l-Z2)) 

* sgnOfl-Z2))  , 

s o dZ/  d p > 0 . □ 
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decreases  with  r increases  with  A and 


Lemma  2.27 . 1 - 


In  Z 


-1 


increases  with  p. 


Proof.  - In  Z = In  Z,  so  the  result  follows  from  Lemma  2.26,  □ 


Suppose  that 


iLlOi  < L ^ Z'V.  A1.  P') 

8(0)-  2 

for  a particular  set  of  parameters  (r',  A',  p’).  Then,  by  the  preceeding 
lemma, 

AliSl  < i _ ln  z'Y^P,). 

g(o)  - 1 2 

for  r < r A>A'  and  p>p'.  Thus,  given  the  ratio  g'(0)/g(0), 
it  is  possible  to  determine  regions  of  the  parameters  ( r.  A,  p)  over 
which  the  hypothesis  of  Lemma  2.24  is  satisfied.  Table  2.1  is  a compila- 
tion of  such  regions  for  various  upper  bounds  on  the  ratio  g'(0)/g(0). 

For  example,  the  first  line  of  the  table  states  that  if  g'(0)/g(0)  < .1, 
then  f(0,0)  > 0 for  r < .3,  A > 1 and  p > .2  (which,  in  view  of 
the  comments  in  15],  probably  hold  when  the  model  is  applicable). 

Recall  that  g( 0)  is  the  total  cost  for  the  first  unit  of  capacity 
added.  Thus,  when  the  expansion  cost  is  concave  g'(0)  is  an  upper  bound 
on  the  marginal  cost  (g(l)  - g(0))  for  the  second  unit  of  capacity  and 


.3 


.3 

.2 

.2 

.2 

.2 

.2 

.2 

.2 


TABLE  2.1;  Regions  where  £(0,0)  > 0 

a lower  bound  on  the  marginal  cost  (over  fixed  cost)  for  the  first  unit 
of  capacity.  Hence,  when  g is  concave, 


Thus,  when  expansion  costs  are  concave,  the  ratio  g'(0)/g(0)  will 
certainly  be  less  than  unity.  In  practice,  the  fixed  (or  "setup")  cost 


I 


\ t 

I 


I 

K 

I ! 


( 


! 


is  typically  much  greater  than  the  marginal  costs,  so  it  is  not  unrea- 
sonable to  expect  g'(0)/g(0)  to  be  relatively  small  when  expansion 
costs  are  concave  (g'(0)/g(0)  < .1  is  not  atypical).  In  summary,  it 
seems  reasonable  to  conclude  that  f will  be  positive  (implying  that  (l£) 
is  satisfied  for  K £ k)  in  most  applications  involving  concave  expansion 
costs. 

Definition  2,7.  For  K £ k,  let  X(K)  denote  the  fixed  points  of  f(*,K) 
over  ( 0,a>)  : 


X(K)  = {X  : X = f(X,K) , X £ (0,»))  . 

If  x(K)  ^ 0,  let  x(K)  = min[X  : X £ X(K)).  If  X(K)  = let  x(K)  = +«. 

Theorem  2.28.  If  holds,  then  |x(K) | < 1,  K £ k. 

Proof.  Let  K £ K.  By  Definition  2.6,  x(K)  ^ (°,  e j = 0.  Also, 

f( e K)  > e and  f'(  • K)  is  concave  over  [e  »).  If  X(K)  = f!>  there 

is  nothing  to  prove.  Hence,  suppose  x(K)  t> , which  implies  x(K)  < +». 

Let  X > x(K)  be  arbitrarily  chosen.  Let  n = (X-x(K))/(x-e^) . Then 
0 < n < 1 and  x(K)  = + ( 1-^)X.  Thus,  by  concavity, 
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= f(x(K),  K)  >tifCeK,K)  + (l.u)  f(X,K) 
> neR  + (l-4)  f(X,K) 

=>X>f(X,K)  . 


Hence^  x(K)  ^ 0 implies  x(K)  = {x(K)). 

Corollary  2.29.  Let  K £ k and  assume  2^  holds. 

(i)  If  X(K)  = 0,  then  f(X,K)  > X,  for  all  X > 0. 

(ii)  If  X(K)  /L  then 


f(X,K)  >X,  X £ (0,  x(K)) 

f(X,K)<X,  X € (x(K),  »)  . 


Proof. 


(i)  and  (ii)  follow  from  Theorem  2.30  and  the  continuity  of 

□ 


Theorem  2.30.  Let  K £ k and  assume  2^  holds. 

( i)  If  f(X,K)  < K,  then  X*(K)  = K. 

(ii)  If  f(U(K)  f K)  > U(K)#  then  X*(K)  = U(K). 

(iii)  If  f(K,K)  > K and  f(U(K),  K)  <U(K),  then 


Proof.  Let  K £ x. 


f 


(i)  If  f(K,K)  < K,  then  X(K)  ^ 0 since  f(e  K)  > 0 and 
f(-,K)  is  continuous.  Hence,  by  Corollary  2.29,  f(X,K)  < X,  X £ (K,co). 

By  Theorem  2.20  this  implies  c*C(X,K)/<3X  > 0,  X £ (K,  U(K)].  Hence 
f(K,K)  < f (X,K) , X £ (K,  U( K)  ] . 

(ii)  If  f(U(K),  K)  > U(K),  then  x(K)  > U(K)  by  Corollary  2.29. 

Hence,  f(X,K)  > X,  X £ (K,  U(K)).  By  Theorem  2.20,  this  implies 
&C(X,K)/SX  < 0,  X £ (K,  U(K) ) . Thus,  f(X,K)  > f(U(K),  K) , X € [K,  U(K)). 

(iii)  If  f(K,K)  > K and  f(U(K),  K)  <U(K),  then  x(K)  £ (K,U(K)) 
since  f(-,K)  is  continuous.  Let  X'  = [x(K)].  Then,  by  Corollary  2 .29 , 
f(X,K)  >X,  X £ [K,  X')  and  f(X,K)  < X,  X £ (X'+l,  U(K)].  Hence,  by 
Theorem  2.20,  dC(X,K)/dX  < 0,  X £ [K,  X')  and  dC(X,K)/&X  >0, 

X £ (X'+l,  U(K)  ].  Thus,  C(X,K)  >C(X',K),  X £ [K,X')  and  C(X,K)  > C(X'  + 1,K), 
X £ (X'+l,  U(K)  ] . Hence,  X*(K)  £ {X',  X'+l).  □ 

The  near-monotonicity  of  X*(*)  is  given  by 

Theorem  2 .31  . If  £ holds,  then  X*(K)  >X*(K')-1  for  all  K,K ' £ x, 

K > K'. 

Proof.  Let  K,K'  £ x with  K > K'.  We  consider  three  cases. 

Case  1:  f(K',K')  < K'.  In  this  case,  X*(K')  = K*  by  Theorem  2.30 
and  X*(K)  > K > K ' = X*(K  ') . 
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Case  2:  f(U(K'),  K :)  >U(K').  In  this  case,  X*(K')  = U(K')  by  Theorem  2.30. 

By  Lemma  2.20,  f(U(K’),  K)  > f(U(K'),  K')  >U(K’).  Hence,  by  Corollary  2.29, 
f(X,K.)  > X,  X £ (0,  U(K'))  and  it  then  follows  from  Theorem  2.30  that 
X*(K)  > U(K’)  = X*(K ')  . 

Case  j:  f(K',K')  >K'  and  f(U(K’),  K)  <U(K').  In  this  case, 

X* ( K ' ) € (ix(K')],  [ x( K ') ] + l}  by  Theorem  2.30.  By  Corollary  2.29, 
f(X,K ')  > X,  X € (0,  x(K ') ) . Hence,  by  Lemma  2.20,  f(X,K)  > X,  X € (O,  x(K')). 
It  then  follows  from  Theorem  2.30  that  X*(K)  > [x(K')j  , which  implies  that 
X*(K)  > X*(  K ’)  - 1,  □ 

With  regard  to  finding  X*(K),  there  are  numerous  methods  of 
computing  fixed  points  for  a function  of  one-variable.  The  technique 
considered  here  is  function  iteration.  It  will  be  demonstrated  that, 
under  Assumption  2.1',  function  iteration  always  determines  X*(K)  when- 
ever condition  holds. 

Definition  2,8.  For  K £ k,  let  d (K)  denote  the  class  of  sequences 
given  by 


XQ  £ [K,  U(K) ] 

X.+1  = max{K,  min(f(X.,K),  U(K))},  i = 0,  1,  2,  ...  . 
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Lerrnia  2.32.  Let  K £ k and  {X^}  £ J(K) . If  g is  concave,  then 

X = lim  X exists  and  X £ [K,  U(K)1.  Furthermore  {X.)  is  monotone 
i — > oo  i 

nondecreasing  (nonincreasing)  if  XQ  < f(X0,K)  (XQ  > f(XQ,K)) . 

Proof.  Let  K£  k and  (X.)  £^(K).  For  each  i,  K<X.<U(K),  by 
Definition  2.8;  hence,  {X±3  <=  [k,  U(K)].  We  first  show  that  (X.)  is 
monotone  by  considering  the  two  cases. 

CaseJ:  XQ  < f(XQ,K).  In  this  case,  Xj  = min{U(K),  f(X0,K)J  > XQ.  If 

x!  = U(K),  then  Lemma  2.22  yields  ffX^K)  = f(U(K),  K)  > f(XQ,K)  > U(K)  = Xp 
since  XQ  £ [K,  U(K)].  Similarly,  if  XL  = f(X  K)  > X then 

f(XpK)  > f(X^,K)  = Xp  by  Lemma  2.22.  In  general,  assume  that  X^  > X.  ^ 
and  f(X.,K)  > X..  Then  X.+1  mintU(K) , f(X.,K) ) > X..  If  X.+1=U(K), 
then  Lemma  2.22  yields  f(X.4,K)  = f(U(K),  K)  > f(X.,K)  > U(K)  = X.+1(K), 
since  X.  £ [K,  U(K)].  Similarly,  if  X.  + 1 = f(X.,K)  >X  then 

f(X£+px)  > f(X^,K)  = by  Lemma  2.22.  Hence,  it  follows  by  induc- 

tion on  i that  (X^j  is  nondecreasing  whenever  X^  < f(X^,K). 

Case  x0  > f(x0,x).  In  this  case,  XL  = max(K,  f(XQ,K)}  < XQ.  If 

X1  = K,  then  Lemma  2.22  yields  f(XpK)  = f(K,K)  < f(XQ,K)  < K = X 

since  XQ  £ [K,  U(K)].  Similarly,  if  *l  = f(XQ,K)  < X0,  then 

f(XpX)  < f(x0,K)  = Xp  by  Lemma  2.22.  In  general,  assume  that  X^  < X^  ^ 

and  f(XpK)  < Xp  Then  X.+1  = max{K,  f(X.,K)).  If  X£  = K,  then 

Lemma  2.22  yields  f(Xi+pK)  = f(K,K)  < f(XpK)  < K = X^,  since 

xt  G fK,  U(K)].  Similarly,  if  X = f(XpK)  < Xp  then 
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by  Lemma  2.22.  Hence,  it  follows  by  induction 


f(Xi+1,K)  < f(XpX)  = Xi+1> 

on  i that  {X^ } is  nonincreasing  whenever  XQ  > f(XQ,K). 

Since  is  a monotone  sequence  and  {X^)  c [k,  U(K) ], 

X = lim  X.  exists  and  X £ [K,  U(K)].  □ 

i -»  oo 

Theorem  2.33.  Let  K € K and  {X^}  £ ^(K).  If  holds,  then 

X*(K)  € [[X],  [X]  + l)  where  X = lim  X.. 

i 7 CO 

Proof.  Let  K £ K and  (X^J  £ J(K) . By  the  previous  lemma,  X exists 
and  X £ [K,  U(K)].  We  consider  three  possible  cases. 

Case  1:  f(K,K)  < K.  By  Corollary  2.29^  this  implies  that  f(X,K)  < X, 

X £ (K,  U(K)].  Hence,  f(XQ,K)  < XQ  and  by  the  previous  lemma,  {X^}  is 
nonincreasing.  Hence,  X ^ = max{K,  f(X^,K)},  i = 0,  1,  2,  ...  . Thus, 
X=  max{K,  f(X,K)}.  If  X > K,  then  f(X,K)  <X  and  max[K,  f(X,K)}  < X, 
a contradiction.  Therefore,  using  Theorem  2.30,  X = K = X*(K)  whenever 
f(K,K)  < K. 

Case  2:  f(U(K),  K)  >U(K).  By  Corollary  2.29,  this  implies  that 
f(X,K)  > X,  X £ [K,  U(K) ) . Hence  f(XQ,K)  > XQ  and  by  the  previous 
lemma,  {X^)  is  nondecreasing.  Hence,  X ^ = min{U(K),  f(X^,K)), 
i = 0,  1,  2,  ...  . Thus,  X = min{U(K),  f(X,K)).  If  X < U(K),  then 
f(X,K)  > X and  min{U(K),  f(X,K)}  > X,  a contradiction.  Therefore, 
using  Theorem  2.30,  X = U(K)  = X*(K)  whenever  f(U(K),K)  > U(K). 
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Case  3 : f(K,K)  >K  and  f(U(K),  K)  < U(K).  In  this  case,  we  consider 

two  subcases  to  show  that  X = x(K). 


Case  j a.  XQ  < ffXQ,K).  In  this  subcase,  xi+ = min{U(K),  f(Xi,K)), 
i = 0,  1,  2,  ...  . Tli us,  X = min{U(K),  f(X,K)}.  But,  f(U(K),K)  < U(K), 
so  X f U(K).  Hence,  X = f(X,K),  which  implies  that  X = x(K). 


Case  3b.  XQ  > f(X0,K).  In  this  subcase,  X = max(K,  f(X  K)J, 

i = 0,  1,  2,  ...  . Thus,  X = max{K,  f(X,K)}.  But  f(K,K)  > K,  so  X ^ K. 
Hence  X = f(X,K),  which  implies  that  X = x(K). 


Tlius,  X = x(K)  whenever  f(K,K)  > K and  ffU(K),  K)  < U(K). 

It  follows  from  Theorem  2.30  that  X*(K)  £ (txfK)],  [x(K)]+lj  = {[X],[X]+1} 
in  this  case.  □ 


The  following  theorem  detai Is  the  convergence  properties  of  the 
function  iteration  approach. 


Theorem  2.34.  Let  K £ k.  Assume  that  >1^  holds  and  that  f(K,K)  > K 

and  f(U(K),K)  <U(K).  Let  {X.j  £ J ( K)  with  x(K)  > XQ  > max{€K,K} . 

If  (X  0 - X )<(X  ,-X)  for  some  m > 0,  then  x(K)  £ (X  , X +1) 
m+2  m+r  v m+1  nr  — ’ v v n>  n 

where : 


ls-1 


n = m+2  + [(In  f,  ) ln(U(K)-K)  + 1] 


and 


l = (X  „ - X .)  (X  . - X )' 
m+2  m+1'  ' m+1  nr 


(2.19) 

(2.20) 


1 


1 
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Proof.  Since  f(X^<K)  > XQ  and  x(  K)  £ (K,  U( K) ) , [X^)  is  strictly 
increasing  and  *1+l  = flX^K) , i > 0.  Since  f(-,K)  is  concave  over 
lX0>“)’  f('»K)  has  nonincreasing  lefthand  derivatives  f'(-K)  over  [X( 
t 4J.  By  concavity,  Leimia  2.22,  and  the  hypothesis 

0 < £:<x„i-K)  5 (f(x„1,K)  - f(V»  (VrV1 

■ (x»2  - Vil  <Vi  - ]iJ'1 

= £ < 1. 

Also,  by  concavity  and  Lemma  2.22, 

x(K)  = f(x(K),K)  < f(Xm+.,K)  + (x(K)  - XJ 

= W + f-(Vi^)  WK)  - 

<X  .1  + ^(x(K)-X  .)  i = 1 2 

— m+i+1  v v ' m+i' 1 > * ••• 

=*X<X>  - x»uiS«*(K)  - *ml),  1.  1,  2,  ... 

- Xm.Ui  ^ W)  - xj,  i - o,  1,  2,  ... 

x(K)  - X < £V(x(K)  - X .) 

' n — v v ' m + 17 


where 


1 


v > (In  i’1)"1  ln(  U(K)-K) 

= |ln  g\~l  ln( U( K) -K) 

= -(In  £)_1  ln(U(K)-K) 

=>  v In  £ < -ln( U(K) -K) 

=»  f < (U(K)-K)-1  . 

Therefore, 

x(K)  - X < (U(K)-K)"1  (x(K)  - X .)  < 1 
11  m+ 1 

=»x(k)  e (xn,  xn+i)  , 

since  U( K)  > x( K)  > X > X > X > K.  n 

n — tth-I  0 — 

With  regard  to  the  condition  on  (X^,  X^^  X^2)  in  the  hypothesis 
of  Theorem  2.34,  the  following  lemma  provides  an  upper  bound  on  m. 

~eJ!Fl2 ..2* 35 • Let  K € <•  As sinne  that  2^  holds  and  that  f(K,K)  > K and 
f(U(K),  K)  <U(K).  Let  {X^  € V(K)  with  x(K)  >max[eK,K).  Then 

(Xmf2  " Xm+1>  < (Xm+1  “ V for 

m=  [(Xi-Xq)-1  (U(K)-K)]  + L . (2.21) 
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Proof.  Since  f(X^>K)  > and  x(K)  £ (K,  U(K)),  {X^J  is  strictly 
increasing  and  X^+^  = f(X^,K)>  i > 0.  Since  f(*,K)  is  concave  over 
fx0j")j  f(’,K)  has  nonincreasing  lefthand  derivatives  f'(.,K)  over 

fXQ>“)  [^1«  ^e  shall  show  that  f'(X^,K)  < lj  the  proof  is  by  contra- 
diction. 

Suppose  fJ(X^K)  > 1.  Then,  since  f(",K)  is  concave  and  (X.) 
is  increasing,  f_'(X.,K)  >1,  i * 0,  1,  ...,  m.  Hence,  by  concavity, 

X.  . - X.  = f (X.  K)  - X. 
l+l  l ' L>  ’ i 


> f(X._1,K)  + f_’(X.,K)  (X.-X.^)  - X. 
>f(Xi_i,K)  + (X.-X.^)  - X. 


= f(Xi-PK)  ' Xi_! 


Xi  * Xi-1  , 


=»X.  - X.  > X,  - xn  , 

i+1  1 — 1 0 * 


^ , 1T1 


i = 1,  2,  ...  , m 


=>  X - Xn  > mfX.-Xj 
m 0 — v 1 O' 


But,  by  (2.21),  the  above  inequality  implies  that  X > X„  + m(X.-X  ) 

m — 0 v 1 O' 

> xn  + (U(K)-K)  >U(K).  This  is  a contradiction,  since  X < lim  X 

’ m , i 

i.  — > oo 

= x(K)  < U( K) . Hence,  we  must  conclude  that  f'(X  K)  < 1 

- m>  ' ' 

Since  f(  ,K)  is  concave  and  {X^  is  increasing, 


1 > f:<V)  > (Vi  - V (f(!wi’K>  - £(V» 

■ <Vi  - v1  - Vi>  • 


Using  Theorems  2.33  and  2.34,  the  Function  Iteration  Algorithm  can 
be  stated  as  follows. 


Algorithm  C , (Function  Iteration,  given  Z^) . 

(1)  If  f(K,K)  < K,  X*(K)  = K (stop). 

(2)  If  f(U(K),K)  >U(K),  X*(K)  = U(K)  (stop). 

(3)  XQ  max  { K,  eR},  «-f(XQ,K),  Xg  *-f(XpK). 

(4)  If  (Xg-X^)  < (XrX0),  go  to  step  (6);  otherwise  continue. 

(5)  XQ  <-Xp  «-X2,  X2  <-f(X2,K);  go  to  step  (4). 

(6)  l <-  (Xg-Xj)  (X^)"1,  X *-Xg,  i <-0,  i*  f-  [ | In  £ | In  ( U(K)  -K)  + 1]. 

(7)  X «-  f(X,K)  . 

(B)  i *-  i+1.  If  i < i',  go  to  step  (7);  otherwise,  continue. 

(9)  X*( K)  e {[X],  [X+l],  [X+2] ) (stop). 

(X0,  Xp  X2)  play  the  roles  of  (X^,  X^p  X^2)  as  given  in  Theorem  2,34. 
By  Lemma  2.35,  the  inequality  of  step  (4)  will  be  satisfied  in  less  than 
[(XrX0)  (U(K)-K)  + 1]  cycles  of  steps  (4)  and  (5);  hence,  the  algorithm 

is  finite.  Step  (9)  follows  from  Theorems  2.33  and  2.34;  X*(K)  can  be 
determined  by  computing  C(*,K)  over  the  three  components. 

By  Theorem  2.31,  additional  efficiency  can  be  realized  in 
Algorithm  whenever  condition  Z holds  and  the  algorithm  is  being  used 
recursively  to  compute  X*( •)  over  k.  Specifically,  for  any 
X'  € K ~ [K,co),  the  initialization  of  XQ  in  step  (2)  can  be  replaced 
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1 

by:  X0<-max{K,  eR,  X*(K')-1}  (since  X*(K ')- 1 < x(K ')  < x( K) ) . 

The  remainder  of  this  section  demonstrates  that  condition  £ wil] 
typically  hold  for  three  important  classes  of  concave  expansion  functions: 
the  concave  power  functions,  the  logarithmic  functions,  and  the  negative 
exponential  functions. 


Example  2.1 . Concave  Power  Functions.  This  example  includes  the  case 
of  g affine, 


g(X)  = 6q  + 1 1(x+i)  > 0,  > 0,  0 < |2  < 1 , 

I 

g'(X)  = Sj  |g(X+l) 


£g-l 


Denote  a = In  Z_1.  Then 

, .11  + Mx+1>  + W 

f(X,K)  = a ln(l  + — i ^ a)  - l 


§ f(X,K)  = 


i1  i2(x+l) 

g2(x+i)g2  - (^(K)+|0)  i^CSg-l) 


6/ 


(^ig)  (X+l)  + a(0(K)  + 50)  618g(X+l)+a|^e2(X+l) 

Taking  one  more  derivative  yields 


!o+l 


2 2| 
sgn(-~  f(x,K))=  sgn((|g-l)  v(X)  - a ||  (X+l)  ^ , (2.22) 

' 3X  * 


where 
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v(X)  , (Ak)  + s0)  + (,2+1)  (X+l,52)  + s^y'2-1 


+ (&K)  + i0)  (^i2)2J  >o  . 

Since  (fi2-l)  < 0;  (2.22)  implies  a2/dX2f(X,K)  < 0 and  hence,  f(.,K) 
is  concave  over  [0,.)  for  all  K £ *.  Hence,  by  Lemma  2.23,  ^ hoids 

whenever  (Tg)  is  satisfied  and  £ holds  whenever  (£^)  is  satisfied 
for  all  K £ k. 

1 

( will  be  satisfied  whenever 
By  Lemma  2.23,  £ will  hold  if 


b s a' 


■ 1 / ^0  + ^l  \ 

K1+-w)  - 


1 > 0 


since  0(0)  > 0 =>  f(0,0)  > b.  Also,  by  Lemma  2.25,  £ will  hold  if 


_ ..  < , a 

8(0)  l0  + £j_  * ?- 


In  particular  (using  Table  2.1),  if  fi^^  + < mlf  then  £ wiU 

hold  whenever  p > .2  and  r < .3. 


— 


Examp le  2.2.  Logarithmic  Functions 


s(X)  = 6-  + 6-  ln(X+l) 


Denote  a = in  Z 


Then 


Taking  one  more  derivative  yields 


'2(X)  = (0(K)  + |Q)  {^(2  lr(X+l)  + 1)  + 0(K)  + | } 


Since  Vl(.)  is  decreasing  and  v£(.)  is  increasing,  ^/^X2  f(X,K)  is 
decreasing.  Furthermore,  v^X)  ->  0 and  v2(X)  -> -ko  as  X -> Hence, 
either 


(a)  Vl(0)  - av2(0)  < 0;  or 

(b)  there  exists  eR  > 0 such  that 


Vj(X)  * av2(X)  / 


> 0. 


< 


X € [o,eK) 


x e (eK,«) 


In  case  (a),  f(*,K)  is  concave  over  [0,oo).  Since  v (0)  > £* 


vx(X)  > av2(X)  =>  ?2  > av2(X)  (X+l) 

> av2(°)  (X+1)  > a^(X+l)  . 

Therefore,  using(2.23),  v^X)  > av2(X)  =»  1 > a(X+l)  f(X,K)  > 1. 

Thus,  either, 

(a')  is  concave  over  [0,«>),-  or 

(b  ) there  exists  > 0 such  that  f(-,K)  is  concave  over 
and  ijf  f(*>K)  > 1 over  [0,  eR). 

In  either  case,  it  follows  that  ^ holds  if  f(o,K)  >0  (and  that  2 
holds  if  f ( 0,0)  >0). 

Hence,  will  hold  if 


I 
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.1  £0  + ^K) 

f(0,K)  = a in(l  + a)  - 1 > 0 . 
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By  Lenina  2.23,  2 will  hold  if 


b s a"  la(  1 + —■  a)  - 1 > 0 , 


since  j6(0)  > 0 =>  f(0,0)  > b.  Also,  by  Lemma  2.25,  2 will  hold  if 


iL!2l  _ 11  < i « 

g(o)  - !0  - 2 


In  particular  (using  Table  2.1),  if  < . 1,  then  2 will  hold 

whenever  p > .2  and  r < .3. 


Example  2.3.  Negative  Exponential  Functions. 


g(X)  = - lL  e 


-l2(x+i) 


t0>h>0’  h>0 


g'(X)  = l2  lj_  e 


-?2(X+1) 


Denote  a = InZ  . Thus 


. -62(X+1) 

-if  *0  ■ ?1  e + ^(K) 
f(X,D  = a 1 In  / 1 + -2  - 


^2  h e 


§ f(X,K)  = 


S2(S0  + 0(K)) 

a(?0  + />(K))  + ^(Sg-a)  e 


-?2(x+i) 
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I 

[ I 

I 


B l 


Thus,  if  < a,  then  d/dX  f(X,K)  is  nonincreasing  and  f(.,K)  is  concave 
over  [0,“).  Assume  Jjg  < a. 

Then  will  hold  whenever 

~i2 

+ PW  \ 

1 > 0 . 


-1  t 5o  - et  c + »m  s 

£(0,K)  . . m(i  ♦ ■■■  — ■ ea ■)  - 


^2  51  6 


The  above  inequality  will  be  satisfied  if 


.1  + **(K) 
a ln(  1 + t~, a)  - 1 > 0 . 


Ml 


By  Lemma  2.23,  Z will  hold  if 


b s a"1  ln(  1 + a)  - 1 > 0 . 
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2U 


*52 


Since  0(0)  >0  and  e < l=>f(0,0)  > b.  Also,  by  Lemma  2.25,  Z 
will  hold  if 


-?2 


&X21  - ^ e.  < i 1 

-T2  - 1 ” 2 


Jo) 


fj0  + ?1  6 


-£2 


Since  0 < e < 1,  the  above  inequality  will  be  satisfied  if 

-1  - - /-  - -------  -1 


?2  ^1  ^0  - 1 ' a/2,  In  Particulsr  (using  Table  4.1),  if  l2  ^ l'1  < . 1, 
then  £ will  hold  whenever  p > .2  and  r < 
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2.5.  Summary 

Given  the  recurrent  stochastic  capacity  expansion  models  pre- 
sented previously  in  [5]  and  [6J,  a Policy  Improvement  algorithm  was 
derived  to  determine  optimal  expansion  sizes  for  general  expansion 
functions.  For  a certain  class  of  expansion  functions,  which  include 
discrete  convex  functions,  it  was  shown  that  each  functional  C(*,K) 
was  unimodal  and  that  the  optimal  expansion  size  X*(K)  + 1 increases 
monotonically  with  K:  an  Interval-Bisection  algorithm  was  given  to 
determine  X*(K)  for  this  case.  The  monotonicity  of  the  optimal  ex- 
pansion size  was  also  demonstrated,  prior  to  integer  truncation,  for 
important  classes  of  concave  expansion  functions  which  are  continuously 
differentiable;  a simple  Function  Iteration  algorithm  was  derived  to 
determine  optimal  expansion  sizes  for  this  case. 


CHAPTER  3 


MODEL  III:  A NON-RECURRENT  GENERALIZATION 

3.1.  Model  III 

Models  I and  II  were  constructed  under  general  assumptions  (a)  - (d) 
listed  In  Section  1.1.  In  order  to  Illustrate  possible  generalization 
for  these  two  models,  assumptions  (b)  - (d)  will  be  simultaneously  relaxed 
for  the  case  of  non-modular  temporary  facilities  (i.e.,  Model  I);  the 
generalized  model  so  obtained  will  be  referred  to  as  Model  III. 

Although  the  generalization  will  be  illustrated  using  non-modular  temporary 
facilities,  similar  generalization  for  the  case  of  modular  temporary 
facilities  (i.e..  Model  II)  is  also  possible;  the  analogies  should  be 
obvious . 

It  should  be  emphasized  that  the  material  given  here  is  preliminary 
in  nature;  its  purpose  is  to  indicate  the  type  of  extensions  that  appear 
to  be  worthy  of  more  detailed  future  investigation. 

The  construction  of  Model  III  will  proceed  under  the  following 
basic  assumptions: 

Assumption  3.1.  All  costs  are  time-stationary. 

Assumption  3.2.  There  exists  a known  integer  $ > 0 which  represents 
the  saturation  population:  an  upper  bound  on  the  level  of  system  demand. 


Assumption  3.3.  When  the  system  demand  £( t)  = i,  arrivals  occur  according 
to  a Poisson  process  with  rate  A^(  i)  and  departures  occur  (independently) 
according  to  a Poisson  process  with  rate  Ag(i)^  = 0 = Ag^). 

Assumption  3.4.  Permanent  facilities  are  always  used  in  preference  to 
temporary  facilities  and  temporary  facility  charges  are  functions  only 
of  the  excess  demand,  the  amount  of  permanent  capacity  available  and  the 
temporary  facilities  usage  limit. 

Assumption  3.1  has  been  implicitly  assumed  in  all  previous  models. 

The  implication  of  this  assumption  is  that  future  advancements  in  technology 
will  offset  inflation  (or  vice-versa). 

Assumption  3.2  differs  from  all  assumptions  used  in  the  construction 
of  previous  models  and  can  be  restated  as  simply  J9(  t)  < jo,  t > 0.  This 
assumption  guarantees  that  the  system  will  never  "explode".  Letting 
Y( t)  denote  the  permanent  capacity  available  at  time  t,  Assumption  3.2 
implies  that  it  suffices  to  bound  Y( t)  < j5.  This  assumption  was 
technically  unnecessary  in  previous  models,  because  the  assumptions  of 
those  models  allowed  characterizations  that  were  independent  of  Y( t) . 

The  physical  interpretation  of  the  parameter  £)  is  obvious:  jg  is  an 

upper  bound  on  the  ultimate  population  of  the  surrounding  service  area  for 
the  system. 

1 
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Assumption  3.3  generalizes  the  demand  process  within  the  Poisson 
framework,  by  allowing  the  arrival  and  departure  rates  to  vary  according 
to  the  system  demand  level.  _If  the  population  of  the  surrounding  service 
area  is  presumed  to  be  constant,  then  this  assumption  can  alternatively 
be  viewed  as  allowing  the  arrival  and  departute  rates  to  vary  according 
to  the  percentage  of  the  surrounding  population  currently  engaging  the 
system's  service.  However,  it  is  important  to  note  that,  given  the 
system  demand  level,  the  arrival  and  departure  rates  are  time- stationary . 

The  restriction  A^(;fl)  = 0 simply  states  that  whenever  the  system 
is  servicing  the  entire  (saturation)  population,  then  the  next  event 
cannot  be  an  arrival.  The  restriction  Ag(0)  = 0 truncates  the  demand 
process;  that  is,  if  the  system  demand  is  zero,  then  the  next  event 
cannot  be  a departure. 

Denote  A(  i)  = A]L(  i)  + A2(  i) , p(  i)  = A1(i)/A(i),  and 

q( i)  = A2(i)/A(i),  0 < i < $.  Then  the  demand  process  can  alternatively 
be  characterized  by  { A( i) , p( i) , q( i) ) as  follows:  if  $( t)  = i,  then 

the  next  event  occurs  according  to  a Poisson  process  with  rate  A( i) 
with  p( i)  the  probability  that  the  next  event  is  an  arrival  and  q( i) 
the  probability  that  the  next  event  is  a departure  [ 3].  Note  that 
P(£)  = q(0)  =0  and  p(0)  = q(jj)  = 1. 

Assumption  3.4  states  that  the  temporary  facility  costs  behave 
according  to  the  non-modular  case.  The  important  point  here  is  that 
these  costs  are  presumed  to  be  independent  of  the  past  history  of  temporary 
facility  utilization.  Given  the  permanent  capacity  available  (Y)  and 
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the  temporary  facility  usage  limit  (K),  the  assumption  presumes  a known 
unique  correspondence  between  facility  costs  and  the  value  of  the  spares 
level  k(  t)  = Y(t)  - £)(  t) . As  demonstrated  previously  in  [5]  and  [6], 
this  assumption  does  not  usually  hold  for  the  case  of  modular  temporary 
facilities. 

Permanent  expansion  costs  will  be  denoted  as  g(X,Y,K)  which 
represents  the  cost  of  adding  X+l  units  of  permanent  capacity  whenever 
there  are  Y units  of  permanent  capacity  available  prior  to  the  expansion 
and  the  temporary  facilities  usage  limit  is  K prior  to  the  expansion; 
as  before,  it  will  be  assumed  that  X > K.  As  in  previous  models, 
expansion  occurs  whenever  $( t)  = Y( t)  + K( t)  (i.e.,  k( t)  = -K( t) ) and 
an  arrival  next  occurs.  When  expansion  occurs  at  time  t,  the  spares 
level  instantaneously  increases  to  k( t+)  = X - K(  t)  and  the  permanent 
capacity  available  instantaneously  increases  to  Y( t+)  = Y( t)  + X + 1. 

Non-expansion  costs  will  be  represented  incrementally.  Fk(K,Y) 
will  denote  the  expected  incremental  (until  the  next  event)  discounted 
costs  while  the  spares  level  has  value  k,  given  that  the  limit  on 
temporary  facilities  is  currently  K and  that  Y units  of  permanent 
capacity  are  currently  available. 

Under  Assumption  3.2,  the  following  bounds  are  nonrestrictive  for 
all  t > 0: 

0 < Y(t)  + K(t)  < fi  , (3.1) 

and 

-K(t)  < k(t)  < Y(t)  . (3.2) 
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3.2.  The  Cost  Transition  Equations 


I 


' 

* 


Let  C^(K,Y)  denote  the  minimum  expected  total  discounted  costs 
starting  from  k(0)  = k,  K(0)  = K and  Y(0)  = Y;  -K  < k < Y,  0 < K < £-Y, 
0 < Y < 5.  For  0 <i  <£>,  denote 

- mki'1  *"d  B(i)  . . (3.3) 

By  conditioning  on  the  time  and  the  type  of  the  next  event,  transition 
equations  for  {C  (K,Y))  can  be  written  in  a manner  analogous  to  that 
used  for  Model  I.  Three  cases  will  be  considered. 

Case  1;  K = Y = 0.  In  this  case,  there  is  only  one  transition  equation 
since  the  only  feasible  spares  level  is  k = 0.  Recall  that  q(0)  = 0; 
hence,  the  first  event  will  be  an  arrival,  triggering  an  expansion.  When 
expansion  occurs,  the  spares  level  will  increase  to  X and  the  permanent 
capacity  will  increase  to  X+l.  Since  is  an  upper  bound  on  the 
permanent  capacity,  it  suffices  to  restrict  the  expansion  size  by 
0 < X < ,ij-l.  The  costs  until  the  first  expansion  are  given  by  FQ(0,0). 
The  first  expansion  cost  will  be  g(X,0,0)  and  the  expected  costs  there- 
after are  (K',  X+l)  (relative  to  the  expansion  time),  where  K' 

denotes  the  temporary  facilities  usage  limit  that  will  be  used  between 
the  first  and  second  expansion  times.  By  (3.1),  it  sufficies  to  restrict 
K'  + X + 1 < fl,  or  0<K'<  jo-X-1.  Whenever  the  first  expansion  time 
occurs,  the  costs  from  that  time  forward  are  given  by  the  sum  g(X,0,0) 

+ Cx  jj(K',  X+l).  We  can  do  no  better  than  minimize  this  sum  over  the 
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feasible  values  for  X and  K'.  Hence^  the  single  transition  equation 


for  this  case  becomes: 


c0(0,0)  = F0(0,0)  + / min{g(X,0,0)  + C (K',  X+l) } p(0)  A(0)  e-A^°) te'rtdt 

0 ^ 


= V°.°)  + ATo)+r  min(g(X,0,0)  + Cx_k(K',  X+l)  } 

= F0(°,°)  + P(°)  min  (g(X,0,0)  + C (K'  X+l)}.  (3.4) 

°<X<£|-1 
0 < K ' < i-X-1 


Case  c.:  0 < Y+K  < j§.  In  this  casef  we  consider  three  subcases. 


Case  2a.  k — Y.  In  this  subcasej  the  first  event  must  be  an  arrival^ 
since  q(0)  = 0.  Conditioning  on  the  time  of  the  first  event  yields 
(recall  that  p(0)  = 1): 


00 

Cy(K,Y)  = Fy(K,Y)  + / Cy_  1(K>Y)  p( 0)  A(0)  e'A(0)  C e"^  dt 


- FV(K’Y>  + 2>ToJ^2i  cy-i<k>1> 


= Fy(K,Y)  + P(0)  cY-1(K,Y)  . 


( 3 . 5a) 


f 


: 

► 


( 


Case  2b:  -K  < k < Y.  In  this  subcase,  the  first  event  is  an  arrival 

with  probability  p(Y-k)  and  a departure  with  probability  q(Y-k) . 

Hence,  conditioning  on  the  time  and  the  type  of  the  first  event  gives 

00 

Ck(K>Y)  = VK>Y)  + £ ( p(  Y-k)  Ck_L(K,Y)  + q(Y-k)  Cfc+1(K,Y)  } 

• A(Y-k)  e-^Y‘k)t  e-rt  dt 

= Fk(K,y)  + a( Y-k)  Ck+1(K,Y)  + 0(Y-k)  Ck  l(K,Y), 

-K  < k < Y . ( 3.5b) 

Case  2c:  k = -K.  In  this  subcase,  the  first  event  is  an  arrival  with 

probability  p(Y+K)  and  a departure  with  probability  q(Y+K).  If  the 
first  event  is  an  arrival,  then  an  expansion  of  size  X+l  (X  > K)  is 
triggered.  If  expansion  occurs,  the  spares  level  increases  to  X and 
the  permanent  capacity  increases  to  Y*X+1.  Proceeding  as  in  Case  1, 
the  minimum  expected  costs  incurred  from  the  expansion  time  forward  are 
given  by 


min  (g(X,K,Y)  + C (K1,  Y+X+l))  . 

K < X < n-Y-1  A_K- 

0 < K’  < J9-Y-X-1 

Hence,  conditioning  on  the  time  and  type  of  the  first  event  yields: 
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C-K<K>Y)  = F_k<K>Y)  +a(K+Y)  C (K,Y) 


+ P(K+Y)  min  {g(X,K, Y)  + CY  „(K  ' , Y+X+l) } 

K < X < a-Y-1  X-K 

0 < K'  < 5-Y-X-l 

(3.5c ) 

The  cost  transition  equations  for  Case  2 are  thus  given  by  (3.5a), 
(3.5b)  and  (3.5c). 


— -S e $ * X+Y  &•  In  this  case,  the  permanent  and  temporary  facility 
capacities  suffice  to  satisfy  the  saturation  population.  Hence,  expansion 
is  never  necessary  in  this  case.  We  shall  again  consider  three  subcases. 

Casejja:  k = Y.  This  subcase  is  identical  to  Case  2a.  Thus 


VK>Y)  = Fy(K,*)  + P(0)  CY1(K,Y)  . 


(3.6a) 


^ase_2b:  -K  < k < Y.  This  subcase  is  identical  to  Case  2b.  Hence 


Ck(K>Y)  = Fk(K>Y)  + a(Y’k)  Ckfi(K,Y)  + P(Y'k)  Ck.i(K,Y)  , 


-K  < k < Y 


(3.6b) 


>c:  k = *K-  ln  this  subcase,  *(0)  = Y-k  = Y+K  = J.  Hence,  since 

p(Ji)  = 0,  the  firsL  event  will  be  a departure.  Thus,  conditioning  on 
the  time  of  the  first  event  yields: 
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C.K(K>Y)  = F_k<K>Y)  +a^)  C-K+1(K>Y)  • ( 3*6c) 

The  cost  transition  equations  for  Case  3 are  thus  given  by 
(3.6a),  (3.6b)  and  (3.6c). 

Using  matrix  notation,  it  is  possible  to  write  a single  unified 
expression  for  the  cost  transition  equations  (3.4),  (3.5)  and  (3.6). 


Definition  3. 

1.  For 

0 < K < £- Y and 

0 < Y < 5,  let 

/P(K+Y) 

min 

K < X < £-Y-l 

(g(X,K,Y)  + CX_K(K',Y+X+1)), 

#k,y)  = < 

0 < K’  < £-Y-X-l 

if  K+Y  < ft 

^0  , 

if  K+Y  = J) 

Using  the  above  definition,  the  cost  transition  equations  (3.4)  - (3.6) 
can  be  rewritten  as: 


C0(°,°)  = F0(°'0)  + #0,0) 

and 

Cy(K,Y)  - 0(0)  C^^Y)  = Fy(K,Y) 

-a(Y-k)  Ck+1(K,Y)  + Ck(K,Y)  - 0(Y-k)  Ckl(K,Y)  = Ffc(K,Y),  -K  < k < Y 

-a(Y+K)  C_r+1(K,Y)  + Cr(K,Y)  = F_R(K,Y)  + #K,Y) 

0 < Y < fl;  0 < K < l- Y,  Y+K  > 0 (3.7) 
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Definition  3.2.  For  m > 0,  let  A(m)  denote  the  square  tridiagonal 


matrix  of  dimension 

mfl  with  nonzero 

elementr 

given 

by 

AiiW 

= 1 , 

i = 0, 

1,  .. 

• y m 

Ai,i+l(nl) 

= -P(i)  , 

O 

II 

1,  .. 

• y in- 

Ai,i-l(m) 

i = 1, 

2,  •• 

. , m 

Let  C(KJY)  denote  the  (K+Y+l)-vector  given  by 


C(K,Y)  = (Cy(K,Y),  Cy_i(K,Y),  ...,  Cq(K,Y),  ...,  C (K,Y))T  . 


Similar ly , let  F(K,Y)  denote  the  (K+Y+l) -vector  given  by 


f(K,V)  - (Fy(K,Y),  Fy1(K)Y),  r0(K,Y),  F_k(K,Y))' 


Using  Definition  3.2  and  the  above  vectors,  the  cost  transition  equations 
(3.7)  become 


A(K+Y)  C(K,Y)  = F(K,Y)  + e(K,Y)  e(K+Y)  , 

0 < K < *-Y,  0 < Y < , (3.3) 

where  e(m)  denotes  the  (mt-l)-unit  vector  with  last  component  unity 
and  all  other  components  zero,  m > 0. 
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Thus,  given  £(K,Y),  the  costs  C(K,Y)  can  be  determined  by 
solving  the  nonsingular  linear  system  given  by  (3.8).  The  next  section 
demonstrates  that  £.(K,Y)  and  C(K,Y)  can  be  computed  recursively. 


3.3.  A Backward  Recursion  for  Computing  the  Minimum  Costs 

In  this  section,  a backward  recursion  is  derived  to  solve  (3.8) 
for  all  feasible  values  of  K and  Y. 


Definition  3.3.  For  0 < Y < £),  define 


7 (Y)  = min  _ C (K’,Y), 
0 < K’  < £-Y 


0 < k < Y 


(3.9) 


Combining  Definitions  3.1  and  3.3  gives 


£(K,Y)  = P(K+Y) 


K < X < 5-Y-  1 
0 < K'  < 5-Y-X-l 


(g(X,K,Y)  + C (K-,Y,.X+l)} 


= f3(K+Y)  min  (g(X,K,Y)  + min 

K < X < fl-Y-1  0 < K'  < ffl-Y-X-1 


min  C (K',Y+X+1)} 


= £(K+Y)  min  (g(X,K,Y)  + Z,  „(Y+X+1) ) , 

X < X < fi-Y-1 

0 < K < fl-Y,  0 < Y < £ . 


(3.10) 


The  recursion  begins  with  the  computation  of  C(0,  fj)  . By 
Definition  3.1,  6(0,  fl)  = 0.  Hence,  from  (3.8), 


1 


(3.11) 


I < 


' . 


1 B i 


-1/ 


c(0,  &)  = A (5)  F(0,  i) 


By  (3.1),  Y = £ implies  K = 0.  Hence,  (3.9)  reduces  to 


\(&)  = ck(°,  = ^.(5)  F(o,  fi)  , o < k < i 


(3.12) 


In  order  to  illustrate  a typical  iteration  of  the  recursion,  let 
Y < & be  arbitrary  and  suppose  that  (Z^Y1),  0 < k < Y * ) are  known 
for  Y+1<Y’  <5.  Then,  for  each  K = 0,  1,  ...,  f)-Y,  g(K,Y)  can  be 
found  using  (3.10)  and  C(K,Y)  can  be  computed  as 


.-1 


C(K,Y)  = A (K+Y)  {F(K, Y)  + £( K,Y)  e(K+Y)J 


(3.13) 


(Zk(Y),  0 < k < Y)  can  then  be  calculated  using  (3.9).  At  this  point, 
the  quantities  {Z^(Y'),  0 < k < Y')  are  now  known  for  Y < Y ' < $ and 
the  next  iteration  of  the  recursion  begins  for  permanent  capacity  level 
Y- 1. 


The  recursion  continues  for  Y = jg-1,  $-2f  ...,  YQ+1,  where  YQ 


denotes  the  initial  permanent  capacity  of  the  system.  Denoting  by  K 


the  initial  temporary  facilities  usage  limit,  G(Kq,  Y ) can  then  be 
computed  using  (3.10).  Finally,  denoting  by  kQ  the  initial  spares  level, 
the  desired  minimum  cost  is  given  by 


.-1 


ck0<K0’  V - \.<VV  <F<K0-V  + s<Ko>V  e<VYo)) 
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The  sequence  of  optimal  actions  resulting  in  the  above  minimum  cost  can 

then  be  determined  by  a forward  lookup  procedure : find  the  sequence  of 

expansion  sizes^  permanent  capacity  levels  and  temporary  facilities 

usage  limits  which  are  the  minimands  in  (3.9)  and  (3.10)  leading  to 

C (K  Y ) . This  is  a standard  procedure  in  dynamic  programming  which 
K0  u u 

is  normally  accomplished  by  storing  all  minimands  along  with  a set  of 
pointers  during  the  backward  recursions  [2  ]. 

R<-  ;all  that  the  system  of  equations  (3.8)  is  tridiagonal.  Given 
the  special  structure  of  a tridiagonal  system^  one  would  normally  not 
compute  the  entire  inverse  matrix  to  solve  the  system  ( see^  for  instance^ 
[5]).  However,  in  the  case  at  hand,  each  inverse  matrix  is  used 
repeatedly;  in  particular^  A ^(m)  is  used  (m  - +■  1)  times  for  each 

m € {Yq+1j  ...j  $}.  llencef  it  may  be  advantageous  in  the  general  case 
to  compute  the  necessary  inverse  matrices  (at  least  in  part)  once  and 
for  allj  prior  to  initiating  the  recursion.  The  following  theorem  pro- 
vides a separate  recursion  for  doing  so. 


5(m) 

= (1  - a(m+l)  p(m)  A^(m))'1  € 1R  , 

(3.14a) 

fc(m) 

= 6(m)  Ct( m+ 1)  A^m);  feT(m)  £lRlfl  , 

(3.14b) 

I3(  m) 

= 5(m)  p(m)  A'*(m)  £®m+1  , 

(3.14c) 

Q(m) 

= A 1(m)  + &(m)  a(m+l)  £(m)  A"^(m)  A^(m)  g ir(  nv+'1)  x( m+1) 

(3.14d) 

Proof: 

By  Definition  3.1, 

A(m)  B(m) 

A(hh-I)  = , 

. c(m)  1 

(3.15) 

where 


B(m)  = (0,  0,  ...,  0,  -P(m))T  £3Rm+1  , (3.15a) 

C(m)  = (0,  0,  ...,  0,  -a(nvfl));  CT(m)  € 1 . (3.15b) 


Also,  A(m)  is  nonsingular  and  so  A_1(m)  exists.  A(nvfl)  is  non- 
singular, so  A \ im-1)  exists;  suppose  that  A^m+l)  is  partitioned 
as  in  (3.14).  Then,  since  A(nn-l)  A \m+l)  = 1,  it  follows  from  (3.14) 
and  (3.15)  that 
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A(m)  G(m)  + B(m)  fe(m)  = I , 


(3.16) 


A(m)  fl(m)  + B(m)  5(m)  = (0,  . ..,  0)T  , 

C(m)  a(m)  + fc(m)  = (0,  . ..,  0)  , 

C(m)  /3(m)  + &(m)  =1 

From  (3.17) , 

A(m)  0(m)  = -5(m)  B(m)  =>#(in)  = -&(m)  A*1(m)  B(m)  . 

Substituting  (3.20)  into  (3.19)  gives 

-S(m)  C(m)  A \m)  B(m)  + 8(m)  = 1 

=*8(m)  = (1  - C(m)  A'l(m)  B(m))"1  . 

Using  (3.15a), 

| 

l 

, A-1(m)  B(m)  = -p(m)  A'*(m)  . 

i 

Hence,  using  (3.15b), 


i 


(3-17) 

(3-18) 


(3.20) 


(3.21) 


(3.22) 


C(tn)  A \ m)  B(m)  = a(m+l)  p(m) 

nrm ' • 


1 


Substituting  the  above  expression  into  (3.21)  gives 

6(m)  = (1  - a(rrH-l)  3(m)  A Vm))  ^ 

mm  ' > 

which  verifies  (3.14a).  Using  (3.16), 

A(m)  Q(m)  = I - B(m)  6(m)  =>  Q( m)  = A_1(m)  - A_1(m)  B(m)  fc(m). 
From  (3.18)^  (3.21)  and  (3.23), 

C(m)  a(m)  + ft(m) 

= C(m)  (A_1(m)  - A_1(m)  B(m)  ft(m))  + ft(m) 

= C(m)  A 1(m)  - (C(m)  A_:L(m)  B(m)  - 1)  ft(m)  = 0 
C(m)  A (m)  = (C(m)  A ^(m)  B(m)  - 1)  &(m) 

^•6(m)  = (C(m)  A-1(m)  B(m)  - l)"1  C(m)  A_1(m) 

= -B(m)  C(m)  A_1(m)  . 

Using  (3.15b), 

C(m)  A \m)  = -a( m+ 1)  ^^(m)  . 


(3.23) 


(3.24) 


(3.25) 
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Substituting  (3.25)  into  (3.24)  yields 


fc(m)  = 5(m)  a(nvfl)  A‘1(m)  , (3.26) 

which  verifies  (3.14b).  Substituting  (3.22)  into  (3.20)  gives 
/9(m)  = 6(m)  P(m)  A'^m)  , 

which  verifies  (3.14c).  Using  (3.22)  and  (3.26)  yields 

A 1(m)  B(m)  fe(m)  = -£(m)  A_1(m)  6(m)  a(nn-l)  A-1(m) 

= -&(m)  a(nH-l)  P(m)  A *(m)  A^1(m) 

Substituting  the  above  expression  into  (3.23)  yields 

G(ni)  ^ A (m)  + &(m)  Ct(nu-l)  (3(m)  A ^(m)  A^*(m)  , 

where  verifies  (3«14d)  and  completes  the  proof.  □ 

In  the  general  case  being  treated  here,  the  minimum  cost  vectors 
C(K,Y)  must  be  computed  for  all  feasible  pairs  (K,Y)  in  order  to  perform 
the  minimization  (3.9).  This  necessitates  a great  deal  of  computation  for 
even  moderate  size  valves  of  £,  as  well  as  substantial  amounts  of  inter- 
mediate storage.  However,  in  many  actual  applications,  it  appears  that 
the  computation  and  storage  requirements  can  be  significantly  reduced. 
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Specifically,  if  the  incremental  expected  operating  costs  for  permanent 
facilities  are  independent  of  the  parameter  K,  then  significant  computa- 
tional economies  accrue.  This  is  demonstrated  in  the  next  section. 


3.4.  A Simplified  Backward  Recursion  for  Semi-Constant  Expected 
Incremental  Operating  Costs 

In  this  section,  it  is  demonstrated  that  significant  computational 
economies  arise  whenever  the  expected  incremental  operating  costs  are  at 
least  partially  independent  of  the  parameter  K. 

Definition  3,4.  The  expected  incremental  operating  costs  will  be  said 
to  be  semi-constant  in  K if  F^(K,Y)  = F^(Y),  0 < k < Y,  0 < K < js-Y, 

0 < Y < £.  The  expected  incremental  operating  costs  will  be  said  to  be 
constant  hi  K if  Fk(K,Y)  = F (Y),  -K  < k < Y,  0 < K < $-Y,  0 < Y < j). 

Obviously,  if  the  expected  incremental  operating  costs  are  constant 
in  K,  then  they  are  also  semi-constant  in  K.  Recall  that  whenever 
k > 0,  the  permanent  capacity  sufficies  to  serve  demand.  Hence,  for 
k > 0,  Fk(K,Y)  will  represent  the  expected  incremental  operating  costs 
for  serving  Y-k  customers  using  permanent  facilities;  it  is  difficult 
to  imagine  practical  situations  where  this  expected  cost  might  depend 
on  the  temporary  facilities  usage  limit  K.  Therefore,  it  appears  that 
the  assumption  of  expected  incremental  operating  costs  which  are  semi- 
constant in  K.  is  extremely  viable  in  practice. 


With  regard  to  the  presumption  of  expected  incremental  operating 
costs  that  are  constant  in  K,  first  note  that  this  presumption  does  not 
rule  out  the  possibility  of  extraordinary  disposal  costs  for  the  temporary 
facilities  just  prior  to  a permanent  expansion.  Since  the  expansion  cost 
function  g is  parameterized  in  K,  these  extraordinary  costs  can  be 
accounted  for  there.  It  is,  however,  possible  to  imagine  situations  wh=>re 
an  increase  in  K will  place  additional  stresses  on  the  entire  temporary 
portion  of  the  system;  in  such  cases,  it  is  hypothesized  that  the  expected 
incremental  operating  costs  Fk(K,Y)  might  increase  with  K for  k < 0. 


Theorem  3.2.  If  the  expected  incremental  costs  are  semi-constant  in  K, 
then  Ck(K,Y)  <Ck(K',Y)  if  and  only  if  Cj(K,Y)  <Cj(K',Y),  0 < j,  k < Y, 

° < K,  K'  < fi-Y,  0 < Y <fi. 

Proof;  Let  0 < Y < 0 < K,  K'  < fi-Y  and  0 < j,  k < Y be  arbitrarily 

chosen.  To  avoid  trivialities,  assume  that  K ^ K ' and  j ^ k.  Without 
loss  of  generality,  assume  k > j. 

Let  Tk(j,  Y)  denote  the  first  time  that  the  spares  level  reaches  j, 
given  that  the  initial  spares  level  is  k and  that  the  initial  permanent 
capacity  level  is  Y.  The  distribution  of  T (j  Y)  is  a function  of 

K. 

the  demand  process  parameters  []\( i) , p( i) , q( i) , i = 0,  1,  ...,  Y-j-1) 
and  is  therefore  independent  of  the  values  K and  K'. 

Let  Gk(j,Y)  denote  the  expected  discounted  costs  over  the  time 
interval  [0,  T (j,Y)].  Since  j > 0,  no  expansion  costs  are  incurred 


74 


j£2 


during  this  time  interval.  Hence,  Gk(j,Y)  is  a function  of  the  demand 
process  parameters  (A(i),  p( i) , q(i),  i = 0,  1,  Y-j-1}  and  the 

expected  incremental  costs  (F^(Y),  £ = j+1,  YJ.  Thus,  G (J,Y)  is 

independent  of  the  values  K and  K'. 

From  the  above  definitions, 

Ck(K'Y)  = GkUY)  + VK>Y)  E[e'rTk(j,Y)]  # 

and 

ck(K',Y)  = Gk(j,Y)  + C.(K’,Y)  E[e“rTk(j>Y>]  . 

Since  k ^ j,  E[e'rX^ J>Y) ] > o.  Hence,  Ck(K,Y)  < Ck(K  ’,Y)  if  and  only 
if  Cj(K,Y)  < CjfK'^Y).  u 

Corollary  3.3.  If  the  expected  incremental  costs  are  semi-constant  in 
K,  then  there  exists  K*(>)  such  that 

Z^Y)  = Ck(K*(Y),  Y)  , 0 < k < Y 

for  0 < Y < & . 

Proof;  For  Y = i,  K*(J5)  = 0 by  (5.12).  For  0 < Y < ®,  define  K*(Y) 
by 

Cy(K*(Y),  Y)  = min  CV(K,Y)  . 

0 < K < J-Y 
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It  then  follows  from  Theorem  3.2  that 

Ck(K*(Y)>  Y)  = min  C (K,Y)  = Z.(Y),  0 < k < Y.  □ 

0 < K < fi-Y 

The  physical  interpretation  of  Corollary  3.3  is  that  there  exists 
a temporary  facility  usage  limit  K*(Y)  that  is  optimal  for  each  Y, 
independent  of  the  starting  spares  level  k > 0. 

From  a computational  point-of-view,  the  above  corollary  is  very 
significant.  The  corollary  states  that  it  is  unnecessary  to  compute  the 
entire  set  of  cost  vectors  (C(K,Y),  0 < K < J-Y)  in  order  to  perform 
the  minimization  (3.9)  determining  tZjjY),  0 < k < Y).  Rather,  it 
suffices  to  compute  only  the  single  elements  (Cj(K,Y),  0<K<iQ-Y) 
for  some  j > 0 in  order  to  determine  0 < k < Y).  In  that 

which  follows,  we  shall  (arbitrarily)  assume  that  j = Y. 

By  (3.8), 

Cy(K,Y)  = A'  L(K+Y)  (F(K,Y)  + e(K,Y)  e(K+Y) ) . (3.27; 

By  Theorem  3.1,  A ^(0)  = [l]  and  for  m > 0 

w • 

A-Vl)  = £0(m)]  , (3.28; 


% (m)  = Ag*(m) 


6(m)  a(mi-l)  p(m)  A'*(m) 


A ( m) 
m.  ' 


#0(m)  = 8(m)  P(m)  AQ^(m)  . 

Hence,  in  order  to  recursively  compute  (Cy(K,Y),  0 < K < ^-Y) 
to  recursively  calculate  {A‘1(K+Y),  A^y  (K+Y),  0 <K  <j9-Y). 
Theorem  3.1,  A'^O)  = [1]  and  for  m > 0, 

= &(m)  5 , 

where 

&(m)  = (1  - a(m+l)  3(m)  A"1(m))"1  , 
fe(m)  = &(  m)  a(m+l)  A'^m)  . 

Combining  (3.28)  and  (3.29),  the  first  and  last  rows  of  A” ^( • ) 
computed  recursively  as  (A~\o)  = [l]): 

5(m)  = (1  - a(m+l)  p(m)  A'^m))"1  , 

mm 

4(m)  = 5(m)  o:(mfl)  A-1(m), 

in. 

A^(nvfl)  = [A^(m),  0]  + A^(m)  (3(m)  [fc(m),  &(m)  ] , 

A^1  = C&(m),  B(m)]  . 


(3.28a) 

(3.28b) 

it  suffices 
By 

(3.29) 

(3.29a) 

(3.29b) 

can  be 


(3.30) 


In  order  to  describe  a typical  iteration  of  the  backward  algoritlan, 
when  the  expected  incremental  costs  are  semi-constant  in  K,  assume  that 
the  first  rows  of  A *"(•)  are  pre-calculated  using  (3.30)  prior  to 
initiating  the  algorithm. 

As  before,  the  recursion  begins  by  solving  the  system 


C(0,  l)  = A_1(fl)  F(0,  5) 


( 3.  31) 


for  Ck(0,  fl)  , 0 < k < jg.  Since  K = 0 is  the  only  possibility  when 

Y = 5,  K*(fi)  = 0 

Let  Y < Jj  be  arbitrary  and  assume  that  (K*(Y),  Y < Y ' < jj) 
are  known.  By  Corollary  3. 3,  Z^Y')  = c (K*(Y'),  Y'),  0 < k < Y 

Y < Y*  < &.  Then  the  iteration  for  Y proceeds  as  follows: 

(1)  K <-  0,  C*  <-  co. 

(2)  fi(K,Y)  <-p(K+Y)  min  (g(X,K,Y)+  Z„  „(Y+X+1) } . 

X < X < £-Y- 1 N 

(3)  Cy(K,Y)  «-A‘\k+Y)  (F(K,Y)  + £(K,Y)  e(K+Y)  ). 

(L)  If  CY(K,Y)  > C*,  go  to  step  (6);  otherwise  continue. 

(5)  c*  <-cy(k,y),  K*(Y)  «-K. 

(6)  K <-K+l.  If  K < $-Y,  go  to  step  (2);  otherwise  continue. 

(7)  Z^Y)  = Cr(K*(Y),  Y),  0 < k < Y. 

In  this  simplified  recursion,  the  partial  cost  vectors  (CY(K*(Y),  Y) , 
. ..,  Cg( K*( Y) , Y))T  must  be  computed  for  each  Y (for  Y = jg  to 
initiate  the  recursion  and  in  step  (7)  of  each  iteration  for  Y < jg) . 
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This  computation  could  be  performed  by  using  a specialized  back- substitution 
scheme  for  tridiagonal  systems  similar  to  that  used  for  Model  I. 


3.5.  Summary 

It  has  been  shown  here  that  many  of  the  assumptions  for  the 
simple  recurrent  models  of  earlier  papers  can  be  relaxed.  For  a given 
saturation  population,  the  revised  Model  III  can  be  solved  by  backward 
recursion  to  yield  an  optimal  sequence  of  expansion  sizes  and  temporary 
facility  usage  limits.  The  price  paid  for  this  generalization  is  that 
of  greater  data  requirements  (arrival  and  departure  rate  estimates  are 
necessary  for  each  demand  level) , increased  computation,  and  the  absence 
of  qualitative  results.  As  shown  in  the  previous  section,  the  backward 
recursion  appears  to  be  computationally  feasible  under  reasonable  assump- 
tions concerning  the  permanent  facility  operating  costs.  It  Is  conjec- 
tured that  the  consideration  of  specific  forms  for  costs  and  demand 
rates  in  the  revised  model  will  prove  to  be  a fruitful  path  of  future 
inquiry. 
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20.  Abstract  \ 

In  two  previous  papers  recurrent  capacity  expansion  strategies  of  the 
form  (X,K)  ''are- consicferedTptfhere  X+l  represents  tha  expansion  size  under- 
taken whenever  excess  demand  (.reaches  the  value  K+l.  The  determination  of 
optimal  expansion  sizes  X*(K)+l'isconsidered  here,  for  fixed  values  K. 

A Policy  Improvement  algorithm  is  derived  to  determine  optimal  expansion 
sizes  for  general  expansion  functions.  For  a certain  class  of  expansion 
functions,  which  includes  discrete  convex  functions,  it  is  shown  that  the 
expected  discounted  costs  are  unimodal  in  the  expansion  size  and  that  the 
optimal  expansion  size  X*(K)+1  increases  monotonically  with  K;  an 
Interval-Bisection  algorithm  is  given  to  determine  X*(K)  for  this  case. 

The  monotonicity  of  the  optimal  expansion  size  is  also  demonstrated,  prior 
to  integer-truncation,  for  important  classes  of  concave  expansion  functions 
which  are  continuiously  differentiable;  a simple  Function  Iteration 
algorithm  is  derived  to  determine  optimal  expansion  sizes  for  this  case. 


It  is  then  shown  that  many  of  the  assumptions  for  the  simple  recurrent 
models  of  earlier  papers  can  be  relaxed.  For  a given  upper  limit  on  the 
ultimate  demand  level,  a revised  model  is  introduced  that  can  be  solved  by 
backward  recursion  to  yield  an  optimal  sequence  of  strategies,  each  of  form 
(X,K) . It  is  shown  that  the  backward  recursion  is  computationally  feasible 
under  reasonable  assumptions  concerning  facility  operating  costs. 


